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Introduction 


Part I of Applied Complex Variables is a largely standard introduction to 
the subject of complex analysis; and I shall not comment on it other than to 
mention that some definitions are not standard. For example, Dettman 
regards a subset to be closed if it is closed as a subset of the Riemann 
sphere (see §1.5 and Theorem 2.1.1). Note that Theorem 2.1.2 would be false 
if he were using closed in the usual sense - consider, for example, f(z) = z 
on the real axis considered as a subset of the plane. However, a closed 
subset of the Riemann sphere is compact. 


In these notes I shall act as a guide to the text and bring to your 
attention points which the reader might overlook; and I shall try to clarify 
the presentation when there seems a possibility of misunderstanding and 
confusion. Not all sections in Chapters 6-10 are in this course and those 
outside it will be clearly indicated. 

To test your understanding there are Self-Assessment Questions (SAQs) and 
Problems. You should try all the SAQs and as many problems as you have time 
for. The latter are sometimes rather hard and so you should not spend too 
much time on any individual problem. Full solutions to SAQs and problems 
will be found at the end of the notes on the relevant chapter. 


In order to pace you through the course we have set four Tutor-Marked 
Assignments (TMAs). These are compulsory in that you cannot pass the course 
without obtaining a reasonable average grade on them. Your three best TMAs 
carry 50% of the total marks for the course, the remaining 50% coming from 
the three-hour examination at the end of the course. Please note that we 
cannot accept any TMAs after the cut-off dates. These cut-off dates are as 
follows. 


TMA 1 (Chapter 6) 

TMA 2 (Chapter 7) 

TMA 3 (Chapter 8, 9) 
TMA 4 (Chapter 10) 


17 April, 1987 
5 June, 1987 

24 July, 1987 

25 Sept., 1987 


Further information about the assessment and the telephone answering set-up 
will be sent to you in the various Stop Presses issued throughout the 
course. 
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It is likely that there are a number of errors in these notes, and that we 
have overlooked some inadequacies in Dettman's presentation. We should be 
grateful if you could inform us of any errors or misprints and of any 
suggested improvements to the Course Notes. 



Chapter 6 
Potential Theory 


§6.1 gives a description of the origins of the mathematics of potential 
theory. The account is not rigorous and the presentation is rather informal 
in places. 


| READ pages 229-232 j 


p.229 The velocity v(x, y, z, t) is the same as the velocity vector v on 

the next page. The physical allusions as to what v(x, y, z, t) means do not 
affect the mathematical status of v(x, y, z, t) which is to be thought of as 
a vector whose value depends upon (x, y, z) in 3-space and time t. 


p.230 In 

Am Am Ax 

—— = —.——. dS = p v dS 
At AV At X 

note that v is the x-component of v and NOT a partial derivative. 

X 


In 


passing from this equation to 
dm 


dt 


= p v.N dS 


it is rather confusing when one tries to think of what 'm' stands for. One 
can think of m as the mass of fluid collected after passing through the 

element as or dS over a period and then ^ is the rate of change of this 

At dt 

mass. Taking account of flow into V and flow out of V through dS one then 
obtains 


dM 

dt 


- JJ ' 

S 


v.N dS. 


p.231 V (nabla) the vector operator, also called 

d _d_ d 
~ d X ' dy r d z 

in cartesian co-ordinates. Hence if v = (v , v , 

~ x y 


v(pv) = ) + + 

~ ~ dx x dy y a z z 


grad, is given by 


v ), then 
z 


The divergence theorem states: 


Note that N points inwards here, but elsewhere in J.W.D. it usually points 
outwards. 


Assuming that + v(pv) is continuous the argument that leads to the 

d t ~ ~ 

continuity equation is the following. If this sum is not everywhere 0 there 
is a point where it is positive or a point where it is negative. By 
continuity there is then a neighbourhood where it is everywhere positive or 
everywhere negative. Taking V to be such a neighbourhood the integral over 
V is positive or negative and not 0. Hence the desired conclusion by 
reductio ad absurdum. 

p.231 A domain D is simply connected if every closed curve in D can be 
continuously deformed within D to a point. If C is a smooth closed curve in 
such a domain then 
f v.T ds = 0 

j <*'W l <"N^ 

c 

is analogous to Cauchy's theorem (see Theorem 3.4.1, p.92). 

To obtain <j> in this case fix (x^, y Q , z Q ) e D and for any other 
(x, y, z) e D let 7 be a smooth path from (x Q , y Q , z Q ) to (x, y, z) and set 
<t> (x, y, z) = J v.T ds. 

7 

One now checks that the integral is independent of 7 for (x, y, z) kept 
fixed and then argues as in the proof of Theorem 3.5.1, p.101 with 
appropriate modifications. 

Self-Assessment Question 

SAQ 1. Which flow is given by the velocity potential 

. v 2 x _ 

<p{x, y) = x + —-^ ? 

x + y 


| READ pages 233-237 j 


p.233 That heat crosses a surface element at a rate 
k vu.N dS 
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is to be taken as a hypothesis of the theory of heat transfer. Vu.N is the 

rate of change of temperature normal to the surface and one is not surprised 
to learn that experiments show that for a given homogeneous substance the 
rate of heat transfer is proportional to this and that the constant of 
proportionality is called the thermal conductivity of the substance. 

Note that heat flows from higher temperatures to lower and so kVu.N dS is 
the heat flow across dS in the direction -N and -kvu.N dS is the heat flow 
across dS in the direction N. Hence if N is inward pointing one sees that 

— = - ff kvu.N dS 
dt JJ ~ ^ 

is correct, i.e. the negative sign jjs necessary. 

p.234, Example 6.1.2 . To say that the curved surface is perfectly insulated 
means that there is no heat transference across it. If u is 'smooth' across 

the boundary then — = ^ =0. A priori one cannot say that u is 'smooth', 

but one argues as if it were and then one checks any 'solution' that one 
obtains. 

p.235. Electrostatics . Mathematically e is a real parameter. If one has 

2 

charges e^, e 2 in place of two equal charges e, then e is replaced by e 1 e 2 * 

K depends on the homogeneous substance in which the phenomena are being 
considered and ^ is called the dielectric constant of the substance. 

A field of force is conservative, by definition, when the work done by the 
force round any closed path is 0. The existence of an electrostatic 
potential for such a field follows in the manner outlined earlier for the 
existence of a velocity potential for irrotational fluid flow. 

p.236, Example 6.1.3 . Mathematically such questions can be considered 
without any special knowledge of the physical context. However, one must be 
able to 'interpret' the latter in mathematical terms. 

2 2 

p.237 v (v u) = 0 is known as the biharmonic equation. 
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Self-Assessment Questions 


SAQ 2. Express v u = 0 in 2 dimensions in polar co-ordinates by considering 
u as the real part of an analytic function and using the Cauchy- 
Riemann equations in polar form. 

SAQ 3. Exercise 6.1, #3. 

SAQ 4. Find the potential in an annulus bounded by circles of radii a < b 
if the inner circle has potential a and the outer potential /? . 


| READ §6.2 | 


p.237. Theorem 6.2.1 The reference to Corollary 3.7.3 in the proof is not 
quite in order. This corollary applies to the real and imaginary parts of 
an analytic function and given a harmonic function it may not be the real 
part of an analytic function in its domain. Consider log|z| in 
{z: 1 < |z| < 2}, for example. Corollary 3.7.3 is applicable to harmonic 
functions in simply connected domains, but does not apply to other domains 
such as an annulus. It is better to apply Theorem 3.7.3 and use that if u 
is harmonic, then u and -u are subharmonic. 

pp.239-240 The function 



2 2 
R - 2rR cos 9 + r 


is known as the Poisson kernel for the disc {z: |z| < R}. It vanishes at 

every point of {|z| = R} except at z = R, where it has a singularity, and 
has the important property of being positive for 0 < r < R. Because of this 
'good' estimates for I , I , I are easily obtained. 

p.239 Note that q(<f>) is continuous on a compact set and so it is uniformly 
continuous. Hence the estimate for |u(r, 9) - g(<9 ) | is uniform in 9, i.e. 
gives e > 0 there is an r(e) (0 < r(e) < 1) such that 
|u(r, 9) - g(0)| < £ (all 9, r(e) < r < 1). 

This observation is necessary in order to establish the first sentence of 
Theorem 6.2.2. 


7 



p.240 The estimates for jl^l, 11 3 1 require 2M in place of M since 

I g (0 ) - g(*)l < |g(^)l + I g ) I < 2 m. 

p.240 For piecewise continuous see Definition 8.1.1 on p.348. 

p.242 Please note that here and subsequently N is the outward pointing 
normal. 


p.243. Lemma 6.2.1 It is better to ignore 
functions v and consider Theorem 6.3.5 and 


this 

f (z) 


derivation of 
R(z - z Q ) 


R - Z 0 Z 


the Green's 


(z, 


x o + 



Z 


x + iy). 


Self-Assessment Questions 

SAQ 5. What conditions on the real numbers a, b, c, d ensure that 

3 2 2 3 

ax + bx y + cxy + dy is harmonic? 


z + 1 

SAQ 6 . By considering arg( -- - - -- .- ) find a function g harmonic in the unit 

disc such that 


g(e 1<? ) = 


1(0 <9 < 7T ) 

-1 ( — IT <9 < 0 ) . 


| READ §6.3 | 

p.245 The essential feature of a Green's function G(x, y, x^, y ) is the 

singular behaviour specified by - —- log|z - z I . The behaviour on C is 

2n 0 

specified so that an application of Green's theorem (Exercises 3.1, #4) 
leads to the desired result. 

pp.251-252 The argument is presented in an unfortunate way. Notation is 
specified subsequent to its use - an undesirable feature in any 
presentation. Perhaps taking one of the two domains to be a disc rather 
than keeping both domains general hinders the understanding. 
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Exercises 6.3, #1 

This seems a rather difficult problem. To show uniqueness we must show that 
if g(x) = 0, then u(z) = 0 in the upper half-plane. Assume then that 

g(x) = o. 

We can extend u to be a harmonic function in the plane by setting 

u(z) = - u(z) for z in the lower half plane. From Green's theorem (Ex 4, 
p.80) we obtain 

2 f (u —) R ae = f (u d -±) R d(9 = f f (UV u + | Vu | ) dxdy. 

J d r J o X. 

0 r=R 0 r=R |z|<R 

Note that vu.N is the derivative of u in the direction of the unit outward 

. . 3 u 

pointing normal N and so in this case Vu.N = —. 


Since u is harmonic v u = 0 and hence 


J < u 

0 r=R 


| |vu| dxdy 
z |<R 


Fixing R n (0 < R Q < 1) we deduce that for R > R n , 


0 


0 


*J |vu| 2 dxdy < fj |vu| 2 dxdy = 2 J (u R ^0 (R oo) . 

'z!<R q 1 z|<R 0 r=R 


Therefore fj |vu| dxdy = 0 and so vu = 0, and consequently u is 

z I 


0 


constant. u is 0 on the x-axis and so u(z) = 0 . 


Self-Assessment Question 

SAQ 7. By considering the definition of a Green's function on p.245 and 
Theorem 6.3.4 prove the result given in Exercises 6.3, #2. 

Problems for s6.1 - §6.3 

PI. Consider a 2-dimensional fluid flow outside a simple curve C in the 

(x, y) - plane. If the velocity v = (v , v ) is such that v - iv is 

x y u y 

analytic outside C, i.e. v - iv = f(z) with z = x + iy and f(z) 

J x y 

analytic, show that the circulation is the same round any simple curve 
surrounding C. 
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P2 . Consider the hypotheses of the preceding problem with C the unit circle 
{|z| = l}. Prove that if the circulation is non-zero then f(z) is of 


the form V a, z with Im a , ^ 0. 
L k -1 

—00 


P3. Find a steady state temperature for a thin metal infinite strip 
{(x, y): -oo < x < oo, 0 < y < l} when the upper edge is perfectly 
insulated and u(-l, 0) = 3, u(l, 0) = 2. How many solutions are there 
for u(x, y)? 

P4. Prove the assertion on p.235 that 


P5. 


P6. 


J 


Ke 


3 


(xdx + ydy + zdz) 


c r 

is independent of the path assuming that the path is smooth. 

If u(z) is harmonic and positive in the unit disc and continuous in the 
closed unit disc show that 

| _ j- iQ 1 *4" IT 

u(0) I - rF < u(re ) < u(0) -.-- - - - (Harnack's inequalities). 

[Hint: first find upper and lower bounds for the Poisson kernel of the 
unit disc.] 

Let D be a domain in the plane bounded by a simple smooth curve C. Let 
g be a real-valued smooth function defined on C. F is the family of 
real-valued functions with smooth derivatives on D and C and boundary 
values g. If 


in£{ JJ 

D 



+ 


2 

u ) dA: u e 

y 


f} 


is attained by U € f show that U is the solution of the Dirichlet 
problem for D with boundary values g. 

(Dirichlet's Principle) 


| READ §6.4 | 


Definition 6.4.1 f(z) would nowadays be called univalent, or perhaps 
schlicht, rather than simple. 


Theorem 6.4.4 The convergence referred to is defined as follows: given any 
compact K c D, f^(z) * f(z) uniformly on K, i.e. given e > 0, there is an n^ 

so that 
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|f n (z) - f(z)| < e (all Z e K, n > n Q ). 

Note that rig depends not only on e , but also on K. 

Theorem 6.4.7 This is one of the most important theorems of complex 
analysis. You should try and get some idea of how the proof goes without 
being too concerned with details. Roughly speaking the steps of the proof 
are as follows. 

(i) Let G be the family of functions univalent (i.e. simple) and having 
moduli bounded by 1 in D. 

(ii) Prove that G is non-empty, i.e. it contains at least one function. 
(This is essential for any subsequent arguments to have substance.) 

(iii) Consider the extremal problem 

M = sup{|g’(z Q )|: g e G} 

with Zq € D fixed. 

(iv) Show that there is an extremal function f € G with |f'(z^)| = M. 

(v) Check that f has the desired properties. 

One of the common reactions to this proof is: Why consider (iii)? If one 
considers a trial and error approach in proving the Riemann Mapping Theorem, 
then one is looking for promising methods that may turn out to be effective. 
If k is analytic then |k'(z Q )| is the local magnification factor of the 

mapping z * k(z) at z = z . The larger |k'(z Q )| is, then the larger the 

magnification at z = z . In the Riemann Mapping Theorem we are trying to 

find a mapping f whose image fills out all of the unit disc, i.e. in some 
global sense we are trying to make the mapping as large as possible. 

However, by fixing a point z^ e D and making the magnification at z^ as 

large as possible subject to (i) one might hope to optimise in a global 
sense, as turns out to be the case. 

Remember that when studying mathematics what we read is often the successful 
outcome of a trial and error approach in the past. The development has 
often been quite chaotic and the orderliness of the way it is presented in 
text-books should be recognised for the sythesis it so often is. 

Self-Assessment Questions 
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SAQ 8. Exercises 6.4, #1. 

SAQ 9. Exercises 6.4, #2. 

SAQ 10. Exercises 6.4, #3. 


| READ §6.5 | 


In dealing with Schwarz-Christoffel transformations note that we know the 
mappings exist by the Riemann Mapping Theorem and Theorem 6.4.8. 

Consequently the existence of the mapping in any particular instance is not 
in question and the problem is to characterise it. This is easier than 
proving its existence ab initio and then determining its properties. It is 
possible to do this without first establishing the Riemann Mapping Theorem 
and Theorem 6.4.8. 

Schwarz - Christoffel transformations are deceptively simple. Given a 
polygon P let z = f(w) be the Schwarz-Christoffel transformation of the 
upper half w-plane onto P. To determine f requires finding the location of 
the points on the real w-axis which map onto the vertices of P. This is 
however an extremely difficult problem even for quite simple polygons P. 

When one deals with Schwarz-Christoffel transformations one can either 
consider the mapping from the upper half plane to the polygon or the other 
way round. It is nearly always better to consider the former mapping and 
not the latter. 

The examples and discussion on pp.260-265 show how one deals with Schwarz- 

dz 

Christoffel transformations, i.e. one considers what must be the form of 

and then finds z = f(w) by integration. The upper half plane has a boundary 
point at w = oo and this may correspond to a point on a side of the polygon 

dz 

or to a vertex of the polygon. If one finds the form of — by considering 

those points on the real w-axis corresponding to vertices, then the correct 
behaviour at w = oo will have been determined thereby without w = oo having to 
be specially considered. 
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Problems for §6.4-§6.5 


I 


P7. Find a univalent mapping of the half-disc (z: |z| < 1, Im z > 0} onto 

the unit disc. 

P 8 . Let f(z) be univalent in {|z| < l} with f(0) = 0. Show that one can 
find g(z) univalent in the unit disc and satisfying 

(g(z )) 2 = f(z 2 ) (|z| < 1 ). 

P9. Let f(z) be analytic and map the annulus {l < |z| < R } univalently 

onto {l < |w| < R^} so that the inner and outer boundaries correspond. 
Show that f(z) = Az for some constant A with |a| =1 and so R = R 2 . 
What happens when {|z| = R } maps to {|w| = l} and {|z| = 1} to 
(|w| = R 2 }? 

P10. Exercises 6.5, #1. 

Pll. Exercises 6.5, #2. 


| READ § 6.6-§ 6.7 | 


3 

P.267 Physically the flows dealt with are considered as taking place in ir 
_ , . . 2 

and when we restrict attention to !R we regard the flow as being independent 
of vertical height. Later in §6.7 when one talks about a 'line source' this 
is to be regarded as a source uniform along a vertical line. 

The formula given for fluid flowing in Example 6.6.1 is in fact the amount 
of fluid per unit height. Note that a similar interpretation must be given 
throughout this section for the results to be thought of as physically 
significant. Flows depending on height are not dealt with. 

Self-Assessment Questions 
SAQ 11. Exercises 6 . 6 , #3. 

SAQ 12. Exercises 6 . 6 , #4. 

SAQ 13. Exercises 6.7, #2. 

§ 6.8 is not part of M828 
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Solutions to Self-Assessment Questions 

1. A uniform stream flowing past a circular cylinder of radius J2. 

2. Let u(z) = Re f(z) and u(z) + iv(z) = f(z). In terms of polars the 
Cauchy-Riemann equations are 

d\l _ dV dV dU 

r Jr ~dl' r ~ 8 r dJ‘ 

Therefore 



8 r 5u x 

8 r r dr ~ I 

and 

so Laplace's 


a , a u 
dr K dr ' i 

3. 

, i 0 v 

^(re ) = 10 

4. 

, i* v (/? 

<Mre ) = — 

5. 

2 

From V </> = 0 


<Mx, y) = 


l ( *i) = J. ( n, 

dr ae 86 dr 


8_ 1 dU 
86 l 86 


£_ ( i £1, 

ae r ae 


log b/a 


3 2 2 3 

= a(x - 3xy ) + d(-3x y + y ) 

3 3 

=aRez +dlmz (z=x+ ly). 

Z 1 

6 . We can define arg( ———) continuously in {|z| < l} so that 


e i<? + 1 

arg(^--> - 

e - 1 


(0 < 6 < n 

{n <6 < 0) . 


Therefore 

, , 2 z + 1 

g(z) = arg (---) 

7T Z -L 

since arg( Z + = Im[log( Z + T )]. 


z - 1 


z - 1 


7. If z = x + iy, z Q = x + iy^, then 

2 2 

log V(x - x ) + (y - y Q ) = log|z - 

2 2 

log \J (x - x ) + (y + y Q ) = log|z - 

and so G(x, y, x Q , y ) is harmonic in (x, 
( V y Q ), (x Q , -y Q ). Also 
G(x, 0, X Q r y Q ) = 0. 


z 0 | = Re log(z - z ) 

z q| = Re log(z - z Q ) 
y) for fixed (x , y Q ) 


except at 
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Assuming that Theorem 6.3.4 is applicable in this case we find that 




u <V V = ■ J 9<x) 


. VG.N dx 


r) 1 7 2 

= - J 9 (x) [- — (- ^ log((x - x 0 ) + (y - y„) ) 


+ log((x - x Q ) 2 + (y + y Q ) 2 ))] dx 

y =0 


An 


J g(x) [- 


2 y 


o 


2 y 


o 


<x - x o )2 + y o 2 (x • x o > 2 + y o 2 


] dx 


y 0 <*> 


= — J 

n •> 


g(x) 


-oo (X - x 0 ) 2 + y 0 2 


dx. 


8 . 


For z , z ^ in {|z| < 1} and z 1 * z^ t 


f(z 1 ) - f(z 2 ) = (z 1 - z 2) (l + a 2 (z i + z 2 * + a 3 (Z l 2 + Z 1 Z 2 + Z 2 2) + * * * ] 
and so 

|f(z 1 )-f(z 2 )|>l z 1 “ z 2 l[!-|a 2 l(l z 1 l + l z 2 l)~l a 3 l(I Z 1 1 2+ I z i z 2 l + I z 2 I 2 )“*•• ] 

> |z 1 - z 21 [1 - 21 a 2 1 - 3 l a 3 l •••] > °- 
Hence if z * z . then f(z ) * f(z ) and so f(z) is univalent in the unit 
disc. 

9. (a) If z e C then the element of length at z along C is |dz|. This 

element of length is transformed by f into an element of length 

along C* given by |f'(z)||dz| at f(z) e C*. Hence 

length of C* = J |f'(z)||dz|. 
c 

(b) If z = x + iy € D, then the element of area at z is dxdy. This is 


transformed by f into an element of area of D* given by |f'(z)| 
dxdy at f(z) e D*. Hence 

2 

area of D* = JJ |f'(z)| dxdy. 


(c) If f(z) = 


k k -1 

a. z (|z| < R), then f'(z) = Y ka z . Consequently 

k 1 ^ 


. i 0,.2 r. . k -1 i (k -1 )6 v - k -1 -i(k-l )6 

(re )| = ka k r e . ^ k^r .e 
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Hence 


area of D* = 


T | f ' (re^ ) | 2 rdrd# 
z I <R 


R 


2 ir go 


J rdr - J <I ka k 

0 0 1 


k-1 i(k-1)0 v - k-1 -i(k-1 )6 


r e 


ka, r e 
k 


) ae 


0 r R V , 2 I I 2 2k " 2 ^ 

= 2n f ) k ja I r . rdr, 

0 1 k 

multiplying up and integrating term by term and using 

r 2 " ik« 1 < k ■ 0) 

e ae = 

0 0 (k *0). 

Therefore 


? ?k 

area of D* = n £ k|a^| R 


> n |a i ! 2 R 2 = n |f'(0)I 2 R 2 . 


2 12 

10.. In place of 3z + z one can consider z + - z . In this case we see 


that in the notation of SAQ 1. 



2 

3 


< 1 and so univalency follows 


from the earlier result. 

11. To say that a flow has a source of strength A at z = z Q means that in 
some neighbourhood of z^ the flow is given by a complex velocity potential 
of the form 

w(z) = A log (z - z Q ) + g(z) (A > 0), 
where g(z) is analytic in {0 < |z - z^| < p } for some p > 0. If C is a 
simple closed curve in this punctured disc that surrounds z = z^ in the 

positive direction then the rate of the flow across C is, apart from fluid 
density, 

t r (- — + g' (z) ) dz = 27 rA. 

l J z - z A 
c 0 

Here we have put J v.N ds into complex form. 

c 


Let f(z) =? be a 1 - 1 conformal transformation of a neighbourhood of z Q , 
i.e. in some neighbourhood f(z) is univalent and so f'(z) * 0 in this 
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neighbourhood. If z - F(f ) and f ^ - f( z g) then F(? ) is 1 - 1 conformal in a 
neighbourhood of ^ and so F 1 (? ) ^ 0 in this neighbourhood. The transformed 
flow at is 

w(? ) = w( F (f ) ) = A log(F(? ) - F(f )) + g(F(?)) 


Since F(? ) is univalent in the neighbourhood of ^ and F (fg) = Zg it follows 
that in some punctured disc {0< |f - ?g| < }r F (? ) ^ z g and so in this 

00 

punctured disc g(F(?)> is analytic. If F(? ) = J a (? - f ) k (| f - ( \ < a) 

k k 0 0 


then 


F (? ) " F (f 0 ) = F(? ) - a Q 


oo a 

■ a i ({ - f o> [1 + \ ^ - V k ~ 1] 


since a^ = F'(f ) ^ 0. Therefore 


log(F(? ) - F(? )) = log (f - f ) + log a. + log[l + 


oo a 

V k , ,k-l 

l r i 


0 " 0 0 ' ”=' ”1 

= log(f - ?g) + function analytic in a neighbourhood 


°f f 0 . 


Hence 


W(? ) = A log(? - ( g) + function analytic in a punctured disc about ? . 

This shows that a source is conformally invariant. 

12. We shall only deal with the essential component of a doublet in a flow, 
i. e. 


w(z) = - 


ic* 

e M 


z - z 


0 


If f = f(z), z = F(f ) is the same as above then 
W(? ) = w(F(f ) ) = - 


ia M 
e M 


F(? ) - F(f 0 ) 

Assuming F(f ) is given as in the earlier solution we see that 
ia 

e M 

w (f ) = “ -—T--r + function analytic at ? . 

«1 If “ f g ; 0 
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Since a = F'(f ) = -- we obtain the result of the question. Setting 

1 U 1 (Zq) 

f * ( z q ) = I f ' ( z q ) I e we see that the direction of the doublet is changed 
from a to 

a + f3 = a + arg f 1 ( z ) . 

13. Considering Theorem 6.2.2 and the result for the integral of the Poisson 
kernel on p.289 we find that if we take z Q = 0, which we can do without loss 

of generality, 


1 r 2 ' , i« . 1 r 2 \ 1 r 2 ’ R 2 - r 2 

_ J u(re )rd* = 57 J [57 J -5-------2 

0 0 0 R - 2r R cos(0 - <p) + r 


is „ 

rd(9 ] u(Re )d <j> 


2n 

= 1 — f r u(Re 1 ^)d^ 

2 , J Q 


Hence 


2n R R 2n 

— J J u (re 1 )rdrd 0 = J rdr. — J u(Re 1 <?!> )d(?!> 
Zn 0 0 0 0 


. u(0) (from Theorem 6.2.2), 


and so 


2rr R 

u( 0 ) = —— J J u(re 1 )rdrd(9 

7T R 0 0 


= ——2 JJ u dxdy (D = {|z| < R}). 


n R D 


Solutions to Problems 

1. Let C^, be two simple closed curves surrounding C in the positive 
direction with one, C say, lying inside the other. Now we can take 


v = (v , v ) and (dx, dy) as T.ds on a curve so that 

~ x y ~ 

v.Tds = v dx + v dy 

~ ~ x y 

= Re[(v - iv )(dx + idy)] 
x y 

= Re[f(z)dz], 

where f(z) is the analytic function v - iv with z = x + iy. By Theorem 

x y 

3.4.5, p.97 it follows that 


J f(z ) dz = J f(z) 


dz 
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and so 

J v.Tds = J v.Tds. 

C x ~ ~ C 2 ~ ~ 

If , C 2 intersect one chooses a curve P that is disjoint from both and 
C 2 and compares circulation round with that round P ; and that round C 2 
with that round r and so one gets the general result. 

2. Since f(z) is analytic in {|z| > 1} it has a representation for such z 
of the form 


f(z) = l 


a k z 


(see §4.6). 


Hence the circulation round {|z| = l} is 

2n 00 k \k() i(9 

Re[J f(z)dz] = Re[J (£ a^2 e )i 2e d0 ] 


z =2 


r . r „k + l r 
= Re[i l a R 2 J 

—oo 


= Re[ 2n ia_^] 


2n 


i(k+1 )0 . 

e d0 ] 


since for any non-zero integer n, 
2tt 

J 

0 


e in ^ d<? = 0. 


Consequently the circulation is non-zero if and only if Im a * 0. 

3. We are looking for a harmonic function u(x, y) in the strip such that 


ay 


u(x, y)| 


= 0 


(—oo < X < oo) 


y=i 


u(-l, 0) = 3, u(l, 0) = 2. 

The upper edge condition suggests looking for u(x, y) which is a function of 
x alone. IF u(x, y) = g(x), then vu = 0 in the strip becomes 

g"(x) - 0 and so g(x) = Ax + B, for constants A and B. Then 
u(-l, 0) = -A + B = 3, 
u(l, 0) = A + B = 2 
and so 


B = j, A = - 
One now sees that 


u(x, y) = - \ + I 
is a solution to the problem. 
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Another possibility is that one looks for u(x, y) of the form g(x) h(y) with 

h'(l) = 0. Consider c real and 

„ r cz-ic. _ r cx+ic(y-l) 

Re[e ] = Re[e J ] 

cx 

= e cos c(y - 1). 

and note that 

cos c(y - 1) = -c.sin c(y - 1) 
dy 

= 0 (y = 1). 

If we choose c = n -, then when y = 0, 

n X tt X 

2 2 n 

e cos c(y - 1) = e cos - = 0 

and so 


7 T X 

U(x, y) = - ^ + | + Ae 2 cos(^(y - 1)) 

solves the problem for any constant A and therefore there are infinitely 
many solutions. 

4. A path C is given parametrically as (x(u), y(u), z(u)) for a < u < b and 
to be smooth means that x(u), y(u), z(u) are nicely differentiable. Then on 
the path 

xdx + ydy + zdz = (x(u)x'(u) + y(u)y'(u) + z(u)z'(u))du 
and so 


Ke 


n ^ v „ 2 r x(u)x'(u) + y(u)y'(u) + z(u)z'(u) 

3 (xyx + ydy + zdz) = Ke J -j - 2 - 2 - 3/2 - 

a (x (u) + y (u) + z (u)) 


du 


= Ke 


o y 2 1/2 

(x z (u) + y (u) + z (u)) 7 a 


Given that one end point of the path (x(a), y(a)f z(a)) is fixed then we see 
that 


J 


Ke 


3 


(xdx + ydy + zdz) 


c r 

depends only on (x(b), y(b), z(b)), the other end of C, and so the result 
follows. 

5. Note first that 


/ ie s 

u(re ) = 


2n 


1 - r 


2 7T 


u(e 1 ^’ )d<?!> 


0 1 - 2 r cos (6 - (p) + r‘ 


and so 
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2 IT 

u(0) = — f u(e 1<f> )6cf> 
2n 0 


i <t> 


Since the Poisson kernel is positive and U(e ) is positive we see that 

2 


2n , . 

1 r / mm 1 - r 

27 J 0<t<2 tt 7 “ " 2 

0 1 - 2r cost + r 


2 tt 

1 r max 

- 2n J 


1 - r 


0<t<2rr „ ■ ^ 2 

1 - 2r cost + r 


irf> i# 

)u(e )d <t> < u(re ) 


) u(e 1</ ’)d <j> 


and so 


i . e. 


1 - r 
(1 + r) 
1 - r 


2n 

1 r . l<t> . , , . 16 . 1 - r 

2 • 2 — J u(e ) d <£ < u(re )< - 

0 (1 - r) 


~ J u(e 1 < ^)d^ 

2 , J Q 


A A 1 - 1 » 

. U(0) < u(re ) < -- r . u ( 0 ) . 


1 + r * ' “ 1 - r 

6 . Let U be a function in F for which the infimum is attained. Suppose 
that v is smooth in D and on C and 0 every where on C. Then for all r? e ir, 
U + r] v e F and so 

JJ (U x + u J) dxd Y ^ JJ t(U + f/v ) 2 + (U + 7] v ) 2 ] dxdy. 


Therefore 


2t) ff (Uv + Uv )dxdy + »/ 2 ff (v 2 v 2 ) dxdy > 0. 

JJ D x x y y J JJ D x y 

If the first integral above is non-zero then by choosing ?? small and of 
appropriate sign we can make the left hand side above negative, which is a 
contradiction. Hence for all such v. 



Consider now the first form of Green's identity in # 4, p.80. With U, v as 
above we see that 
2 

M (W u + vv.vu)dxdy = J vvU.Nds. 

D c 

Since v = 0 on C and JJ vv.vu dxdy = 0 we see that 

D 

2 

JJ W U dxdy = 0. 

D 
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2 

Suppose now that v U is continuous in D and at some point (x^, y^) € D is 
non-zero. Then for some e >0, d > 0, {|z — z | < 2d} e D, where 


= x 


+ iy , and V U > e in this closed disc or V U < -e in the disc. 


Define v by 

1 (|z - z Q | < d) 

v(z) = (r - d ) 3 ( 2 d - r ) 3 (d < |z - z Q | < 2 d, |z - z Q | = r) 
0 ( | z - z Q | > 2 d). 


This particular v satisfies the requirements of smoothness and being 0 on C 
and so 


2 2 
0 = |JJ v V U dxdy| = I J v V U dxdy 

D z-Zg|<2d 


> £ 


[I 


dxdy 


z-z Q |<d 


= e . n d > 0 . 

2 

This is a contradiction and so V U = 0 in D, i.e u is harmonic in D and 
therefore solves the [Dirichlet] problem for D with boundary values g. 

7. Consider §2.4, p.43 on Mobius (linear fractional) mappings. In this 
problem the boundary of the domain consists of two circular arcs if we think 
of a straight line as a circle through z = oo. The picture is clearer if we 
consider the Riemann sphere rather than the complex plane. If we employ the 
1 z 

mapping w = ---, then in the w-plane the image will be bounded by arcs of 

circles through z = oo, i.e. straight lines and intersecting at z = 0 . 
Furthermore the angle between the arcs of the boundary at z = -1 is - and so 


the image arcs or segments will also be at right-angles. 


Clearly the segment between -1, 1 of the real z-axis maps onto the positive 
real w-axis. By conformality, or by direct verification, the semi-circular 
maps onto the positive imaginary w-axis and so the half-disc in the z plane 
maps onto the first quadrant of the w-plane. 


2 

The mapping W = w maps the first quadrant in the w-plane onto the upper 

W - i . 

half W-plane and ? = - + v - maps the upper half W-plane onto the unit f-disc 
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about ? =0. Combining these mappings one finds that a mapping of the kind 
asked for is 


/'N 



\ 


<r 


W - i 
W + i 


2 

w - 1 


2 

w + 1 


( 


( 


1 + z 2 
1 - z' 

1 + z 2 

1 - 7 ? 


+ 


i 


i 


2 2 

( 1 +z) -i(l-z) 

2 2 * 

( 1 +z) + i (1 - z) 


8 . If f(0) = 0/ then f(z) = z F(z), where F(z) is analytic in {|z| < l}. 
Since f(z) is univalent in the unit disc F(z) has no zero in {|z| < l}. 

F(z) cannot be 0 for z * 0 and if F(0) = 0, then f'(0) = 0 and then f(z) 
takes values at least twice in a neighbourhood of z = 0. (See Theorem 
6.4.1). Hence F(z) has an analytic logarithm in {|z| < l}, i.e. there is a 
function h(z) analytic in the unit disc such that 


F(z) = e h(z) (|z| < 1). 

Consequently, F(z) has an analytic square root 

„ , „ .h(z ) . 

G ( z ) = exp [—--] 


since 


(G(z)) 2 = e h(z) = F(z). 

2 

Consider now g(z) = z G(z ). g(z) is analytic in {|z| < l} and 

2 2 22 22 _ 2 

(g(z ) ) = z (G(z )) = z F(z ) = f(z ). 

2 2 2 

Also G(z ) is even, i.e. G(z ) = G((-z) ), and g(z) is odd, i.e. 

g(-z) = -g(z). Suppose g(z) were not univalent in the unit disc. Then for 


some z , z 2 in {|z| < l}, z * z 2 , g(z ) = g(z ) and so f(z 1 ) = f(z 2 ). As 


f(z) is univalent we must have z^ ~ Z 2 ' Hence z i = Z 2 or z i = _Z 2 ’ 

Assuming as we are that z, * z„, we see that z, ^ 0 and z„ = -z,. Hence, 

12 1 2 1 

since g is odd, 

g(z x ) = g(z 2 ) = g(-z l ) = -g(z 1 ) 

and so g(z^) = 0, z^ * 0. However, from the above one sees that the only 

point in the unit disc where g(z) = 0 is z = 0. From this contradiction the 
univalency of g(z) follows. 

9. Considering Corollary 4.5.2 and Exercise 4.5, # 5 one arrives at the 
following result. If f(z) is a univalent map of {l < |z| < } onto 

{l < (w| < R 2 ) with the inner and outer boundaries corresponding to one 

2 

another one can continue f(z) into {R^ < |z| c } by setting 
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? 1 2 

f(z) = 3 - (? = r R, < |z| < R-, )• 

f(0 2 

In this way one can continue f(z) to give a univalent mapping of {|z| > 0} 
onto {|w| > 0 } with f(z) being bounded near z = 0 ; in fact, f(z) -» 0 (z -* 

0). Hence f(z) has a removable singularity at z = 0 and if we assume it has 
been removed then f(z) is entire (Definition 3.6.1, p.115) and univalent. 


By Picard's Theorem (Theorem 4.6.5, p.167 with z^ = oo) we see that 

f(z) = Az + B and since f(0) = 0 we must have that f(z) = Az. The unit 
circle is mapped onto the unit circle and hence |A| = 1. Consequently 

R 2 R l* 

10. (Consider also Example 6.5.5, p.262) 

From the point of view of conformal mapping the strip has two boundary 
points at oo according as we approach oo from the right or from the left. On 
the Riemann sphere there is one point at oo, but in conformal mapping a point 
on the boundary in the set theoretic sense has to be considered together 
with how we approach the point from within the domain in question. 


Let z = z(w) map {im w > 0} onto the strip so that z(l) = +<*> and z(-l) = -<*> 
where too, -oo is oo approached from the right and left respectively. By the 
Riemann Mapping Theorem and Theorem 6.4.8 we know there is such a mapping 
f (w) with $ , on the real axis corresponding to -oo, +oo, i.e. f U x ) = ~°°> 

f ($ ) = +oo. Then for suitable real constants A > 0, B we have that 


z(w) = f (Aw + B). 

From the form of the Schwarz-Christoffel transformation 

dz c 

dw 2 

1 - w 

for some constant C. Hence 

c - 1 + w 

z = - log -- + d. 

2 3 1 - w 

Considering that as w traverses a small semi-circle from 1 - e to 1 + e in 
{im w > 0} Im z increases from 0 to n we see that c = 2. Since the strip is 
invariant under lateral translations, apart from being real d is arbitrary. 
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Taking d = 0, then 


z 1 + w e Z - 1 

e = --, w = - 


w 


e Z + 1 


Also 


<r = 


w + i 

is a univalent mapping of {im w > 0} onto {|f | < l}. 

* - 4' 11 - 1(e z + 1} - f < Z > 

(e - 1) + i(e +1) 
maps the strip onto {|? | < l). 


Consequently 


To say that g is continuous on the boundary of the strip means that g(x), 
g(x + in-) are continuous for -<x> < x < oo and that 
lim v lim 

x -*» g < x > = g(* + = g (+«»), 

lira , , lim 

X+-00 9(X) = x-^-oo 9(X + 

where g(+<x>), g(-oo) are finite. If z = F(? ) is the inverse of ? = f(z) then 

g(z) = g(F(f )) = G(? ) is the function on {|?| = 1} corresponding to g(z) on 
the boundary of the strip. Let u(? ) be the solution of the Dirichlet 
problem for {|?| < 1} with boundary values G(f). When u(f(z)) = U(z) is the 
solution of the Dirichlet problem for the strip with boundary values g(z). 
11. One can adapt the result of Example 6.5.4 or one can argue as follows. 
The (half) strip has 3 boundary points — 0 , n , oo. From the Schwarz— 
Christoffel theory one know that {im w > 0} be mapped onto the half—strip so 
that a, b , c on the real w-axis correspond to 0, ^ oo for some a, b, c. By 
supplementing with a Mobius transformation one can assume that a, b, c are 
0 r If oo. Hence the transformation required is of the form 
A 


dz 

dw 


W 1/2 (l 


w) 


1 / 2 * 


The segment [0, 1] of the real w-axis maps onto the segment [0, *] of the 

d z 

real z-axis and hence is real for 0 < w < 1 and so if we take 


w 1 / 2 (1 


1/2 


w) to be positive for 0 < w < 1, then we see that A is real 
(and positive). Integrating 

A A 


dz 

dw 


(w(l - w)) 


1/2 n 1 ! 2 , 12 1/2 

((j) " ( w ~ j) ) 
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we obtain 


„ . —lw - 1/2 v 

z = A sm - ) + B * 


Now 


. -lw - 1/2 

sin ( 172— ) = 


(w = 0) 


(w = 1) 


and so 


or 


. -1 W - 1/2 n 
z - sm + 2 


n w - 1/2 

Sln(z - 2 “ T75- 

1 l 

- - w = - COS z. 

2 2 


Hence 


1 . 2z 

w = -(1 - cos z) = sin -. 
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Chapter 7 

Ordinary Differential Equation 


| READ §7.1 | 


p.281 


It is logically preferable, I think , 

oo n oo n 


l ~r cos n(0 - <f>) 
n=0 R 


£ 

Re[ Y — (cos n( 
n=0 R n 


n=0 


to consider for 0 < r < R 
- cf>) + i sinn(0 -<£))] 


t 


p -f z 

The rather abrupt introduction of -———- is to some extent 'putting the cart 
before the horse'. 

p.282 Note that x = p sin <f> cos 9 , y = p sin <fi sin 9 , z = p cos <p and so we 

a 2 u 

can find the derivatives of u with respect to p , <f>, 9 in terms of —, —- 

d x 

and so on. Hence Laplace's equation in spherical co-ordinates can be 
verified. 


Self-Assessment Questions 
SAQ 1. Exercises 7.1, #1. 

SAQ 2. Exercises 7.1, #2. 


SAQ 3. Verify that the Legendre polynomial 

n 

w(z) = P (z) = —-- [(z - l) n ] 

n ^n , _ n 

2 n! dz 

satisfies the Legendre equation 


(1 - z )W" - 2zw' + n(n + 1)W = 0 

r 2 n 

[Hint: Put v = (z - 1) and show that 


(1 - z )v" + 2(n - l)zv' + 2nv = 0 and differentiate n times.] 


Preliminary discussion to §7.2 

Given an ordinary differential equation and an initial condition it may not 
be clear if there is a solution to the problem or if there are several 
distinct ones. To cope with such problems we require existence and 



uniqueness results. The simplest equation is of order 1 and it is with such 
equations that Theorem 7.2.1 deals. 


An equation of order 1 is of the form 
d v 

0(x, Y) dx + ^ (X ' Y) = 0 
or 


^ = f(x, y) (f(x, y) 

One is assuming 0(x, y) ^ 
equation for Theorem 7.2.1 

x(t) = — = u(t, x). 


0(x, y) 

0 (x, y) ' 

0 and this seems a reasonable assumption, 
is this one in a different notation, i.e. 


The 


A very simple equation of this form is when u(t, x) = tx and then, 

t 2 /2 . . 

separating the variables, x(t) = Ae for some constant A. An initial 


condition x(t Q ) = x Q 


will determine A uniquely from x^ = Ae 


t o 2/2 


What seems a simple equation is obtained by taking u(t, x) = t + x. 
However, obtaining an explicit solution in this case is not at all trivial. 


If one assumes a solution of the form x(t) 



one can determine the a, 

k 


and then the solution is 

2 

x(t) = -2(1 + t + fj) + (2 + A)e\ 

for some constant A. An initial condition will determine A and so give a 
unique solution. 


2 2 

Taking u(t, x) = t + x leads to an equation that appears not to have any 
nice explicit solutions. However, proceeding as with the preceding equation 
one can show with some difficulty that a solution does exist. One sees that 
what seems a very nice function for u(t, x) can lead to an equation which is 
not easy to handle. The result of Theorem 7.2.1 gives conditions on u(t, x) 
that enable one to say, at any rate, that there is a solution and for given 
initial condition it is unique. It is necessary to impose some conditions 
on u(t, x) so that the problem becomes manageable and those of Theorem 7.2.1 
are weak enough so that for most applications the result is adequate. 
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Existence and uniqueness theorems are of intrinisic interest to the pure 
mathematician, but even for applications they can be useful. It may be 
possible to obtain information about some physical problem if one knows a 
unique solution exists, from the equation itself perhaps, without having any 
explicit form of this solution available. 

In §7.2 one finds that having coped with the first order equation one can 
extend the results to higher orders. Therefore Theorem 7.2.1 is of 
importance that extends beyond its immediate context. 


| READ pp.284-287 | 


Theorem 7.2.1 The proof proceeds by solving the integral equation 

t 

(*) X(t ) = X + f U[r , X(r ) ] dr 

and showing that this equation is equivalent to the differential equation 
with initial condition. The integral equation has two advantages over the 
initial form of the problem: 

(i) the initial condition x(t^) = x^ is built into (*); 

(ii) it is much easier to estimate the effect of integration on a function 
than the effect of differentiation. 

The approximation argument of the proof is quite 'natural' and you should 
avoid being confused by the detail. Note two essential requirements: 

(i) at each stage one must check that x (t) satisfies the condition at 

K 

the top of p. 285 so that x (t) can be properly defined; 

Kt 1 

(ii) one must check that x (t) converges to a limit function with the 

K 

desired properties. 

As far as (ii) is concerned, because the argument proceeds by successive 
approximations it is appropriate to think of the limit function as given by 

OO 

X 0 + I < x n (t) - x n-l (t >>- 
n=l 
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P» 287 The formula preceding Theorem 7.2.2 should begin z = x + iy and not 


z = x + iy. 


| READ pp.287-296 | 


Pt 288 The definitions and properties of functions analytic in two variables 
and associated double power series are dealt with earlier in §4.7. 

P•289 The integral equation at the top of the page is arrived at on the 
supposition that the differential equation of Theorem 7.2.3 can be solved. 

It is, therefore, an alternative form of the problem. 

&J-2.90 Note carefully the use of majorisation in the 'Second Proof' of 
Theorem 7.2.3. 

Self-Assessment Questions 
SAQ 4. Solve the equation 

x(t) = t 2 + x 

with initial condition x(0) = 1. 

SAQ 5. Show that for any prescribed initial condition, 

. 2 2 
x(t) = t + x 

has an entire function solution. 


| READ §7.3 | 


Definition 7.3.1 Note the linear dependence for functions requires 
c ]L w i (z) + c 2 w 2 (z) being 0 for all z in {|z| < R}. It is essential that 

later w^, w 2 are thought of as functions and not as special values. 

Theorem 7.3.1 I find the proof of this theorem rather unconvincing and 
perhaps the following outline is more persuasive. If w^, w 2 are linearly 

dependent, then for each z^ in {|z| < R] the pair of equations at the top of 
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p.298 have a non-trivial solution and so the determinant of the system VJ(z Q ) 
is 0. Hence W(z) = 0 in {|z| < R}. 


Suppose next that W(z) = 0 and that w 1 ( z q) * O' f° r some z q i n U z l < R)• 
There is a neighbourhood of z^ in which w^(z) * 0 and for such z 


A W (Z) 

—- ) = 

dz v w^(z) 


W(z) 

2 , 
w (z) 


= 0 . 


In this neighbourhood, therefore, w 2 ( z ) = c.( z ) for some constant c. By 
analytic continuation this equation holds throughout {|z| < R} and so w^, 
are linearly dependent. If there is no such z , then w^z) = 0 in {|z| < R} 
and in this case linear dependent is trivial. 


Problems for §7.1-57.3 
PI. Exercise 7.1, #3. 

P2. Exercises 7.2, #1, excluding final part of question. 
P3. Exercises 7,3, #1. 

P4. Exercises 7.3, #2. 

P5. Exercises 7.3, #5. 


| READ §7.4 | 


Theorem 7.4,1 It may not be clear in the proof where the condition 

Re[f(0)] > 1 is used, since this is not pointed out. From the recursive 

relations for c, and C, one observes that each coefficient will be properly 
k k 

defined if and only if k + a^ 0 for k > 0. As a^ = f(0), the condition 
Re[f(0)] > 1 ensures that k + a^ ^ 0. 

The condition Re[f(0)] > 1 is considered because of the discussion preceding 
the statement that Re[f(0)] ^ -k for any non-negative integer k. 

p,309 The Green's function considered here is rather different from those 
considered in the preceding chapter. However, in each case G(z, ? ) is a 
kernel which if integrated appropriately together with a given function 
gives a function of the kind sought. In the earlier chapter the function 
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sought was a harmonic function with prescribed boundary behaviour and in 
this case a particular solution of the non-homogeneous differential 
equation. 

Problems for §7.4 

P.6, 7, 8, 9, 10, 11 are Exercises 7.4, #1, 2, 3, 4, 5, 6. 


| READ §7.5 | 


In this section factors like z^ occur where v may not be an integer and so 

V 

z may have different determinations. In such instances one should think of 
v 

z being positive for z small and positive. Other values are obtained by 
analytic continuation. 


00 00 

p.312 In the formula for J (z), £ should be £ . 

U n=l n=0 

P.312 The discussion about the linear independence of J (z), J (z) when u 

v —v 

is half an odd positive integer is prompted by the fact that in this case 
v - (~u) = 2v , a positive integer. However, it may happen that the series 
method for solving a differential equation will yield linearly independent 
solutions even though the roots of the indicial equation differ by a 
positive integer and this is the case here. 

Note also that from the form of Bessel's equation if J (z) is a solution 

v 

then so is J (z). When v is half an odd positive integer then J (z), 

-v V 

J (z) behave differently near z = 0; one is bounded and the other is 
-v 

unbounded and so they must be linearly independent. 

Problems for §7.5 

P.12, 13, 14, 15, 16, 17 are Exercises 7.5, #1, 2, 3, 4, 5, 6. [In #5, y' 
should be y.J 


| READ §7.6 | 
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pp.318, 319 The natural approach would be to find the form of P^(z) by 
considering a solution of the form 
p(z) = a Q + a ± z + - - a n z n 


and finding a recursive relation that specified each of a , a , 
of the preceding. However, we can also consider 


in terms 


q(z) = 2 n ,a n + - - -) . z n (b 0 + ^ + - - A) 

z 

and find each of b , b^,... in terms of the preceding. For the Legendre 


polynomials, if one 'knows the answer' then one realises that is is better, 
i.e. simpler, to use q(z) rather than p(z). 


Problems for $7.6 

P.18, 19, 20, 21, 22 are Exercises 7.6, #1, 2, 4, 7, 8. 


j READ pp.325-333, excluding Theorem 7.7.4 | 


Sturm-Liouville theory has a number of important applications, but a 
rigorous account of the theory is rather difficult. Dettman's presentation 
has some drawbacks and I hope the following commentary will help you to cope 
with these. 

p.327 et seq. The properties of the Wronskians that arise in this context 
are not adequately dealt with. If u, v are solutions of 
y" + [A - q(x)]y = 0 

then 

d 

W(u, v) = uv" - u"v = —u.( A - q) v + (A - q)u.v = 0 

and so W(u, v) is constant, i.e. independent of x. This is observed at the 
top of p.329. Hence one sees that the result about the Wronskian and linear 
independence considered earlier is still true here. 

Theorem 7.7.3 As far as proving that <f>, </’ as defined in Theorem 7.7.2 
satisfy the given integral equations one could do this by direct 
verifications and so avoid the iteration procedure. This procedure does 
seem necessary in order to prove that W(A) is an entire function of A. 
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p.331 Near the top of the page the claim that </>(x, A), ip(x, A) being entire 
implies that W(A) is entire is not quite correct since W(A) also involves 
<p'(x, A), ip'{x, A). From the differential equation one sees that 4>", ip" are 
entire in A since <p , ip are and hence it follows that <p ' , ip' are. 

p.332 The reference to 'the Fourier series of f(x)' near the bottom of the 
page is a bit misleading. The series arising in Sturm-Liouville theory are 
a generalisation of Fourier series involving sines and cosines. 

Theorem 7.7.4 In the proof of the theorem note that cos ? (x - a) means 
cos[f (x - a)] and not (x - a).cos ? . 

You may not find the formulas at the beginning of the proof as easy to 
establish as the author assumes. The first of these can be obtained as 
follows. 

Let <p (t, A) be the solution of <p " + (A - q) <p = 0 with (p { a) = sin a, 

2 2 
<P ' ( a) = - cos a and A = ? . Then, with <p (t) = <p (t, ? ), 

x x 

J sin (x - t)}q(t) 0(t)dt = J sin {?(x - t)}[0"(t) + ? 2 0(t)]dt 
a a 

and 

x xx 

J sin{? (x - t )}<£"( t) dt = sin{? (x - t)}^'(t)| + ? J cos{<r (x - t) }0 ' (t) dt 
a a a 

x 

= -sin{? (x - a)}<£'(a) + ? cos{? (x - t)}^(t)| 

a 

2 X 

- f J sin{? (x - t) }0 (t)dt 
a 

= cosa . sin{? (x-a) } + ?|>(x) _ sine*. cos{? (x-a) } ] 

2 x 

- ? J sin{? (x - t) }0 (t)dt. 
a 

2 

Using this result in the initial equation gives the formula for <p{x, ? ). 

p.335 Figure 7.7.1 is relevant to the particular case a = 0, b = 1. For 

given parameters a, b one has to replace (n + ^)tt by (n + ^)-— . 

2 2 b - a 
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P.336 The verification that the formula at the bottom of the page follows 
from the preceding is quite lengthy. One changes the order of integration 

x b 

and uses the calculus of residues on each of the integrands in J and J . 

a x 

If these two integrals are then evaluated their contributions fit together 
nicely to give the final formula. 

P.337 The concluding sentence of the proof of Theorem 7.7.4 may leave you 
confused. The discussion preceding the statement of Theorem 7.7.4 shows 
,, , lim 1 r 

that 2~n~i J Y( x ' A )dA is expressible as a series with respect to the 

°o 7T i c 

n 

eigenfunctions. In the proof of Theorem 7.7.4 it is shown that this limit 
is also expressible as Fourier series. Using results about the convergence 
of Fourier series from §8 one establishes the result of the theorem. 


P.338 The boundary conditions in Corollary 7.7.4 are 
cos a . y(a) + sin a.y'(a) = 0, 
cos 0 . y(b) + sin /J.y'(b) = 0, 

as in (a) on p.326. Noting on p.332 that ?>(x, A ) = A <P(x, A ) and on 

K K K 


k 1/2 

p.333 that (x) = [ — ■ - —-] <P(x, A, ), then, since <p { a, A, ) = sin a, 

K w (A ) k k 

k 

<£'(a , A ^) = - cos a and ip (b, A ^) = sin 0, ip ' ( b, A ^) = - cos 0, it follows 

that the integrated term on line 3 from the top is 0 at both upper and lower 
limits. 


Problems for §7.7 

P.23, 24 are Exercise 7.7, #1, 2. 


| READ §7.8 


If one assumes that G has properties 1, 2, 3, 4 on p.340, then one can prove 
directly that the solution of (py')' - qy = -f satisfying the boundary 
conditions 

cos a . y(a) + sin a.y'(a) = 0, 
cos 0 . y(b) + sin /?.y'(b) = 0, 


is given by 
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b 

y(x) = J G(x, t) f(t)dt. 
a 

Consider G = G(x, t) with x fixed in (a, b) as a function of t and 

G' = G'(t) = j^-G(x, t). Then 

b x b 

J G(py')'dt = (J + J ) G(py')'dt 
a ax 

x x 

and, thinking of f = im f with <r < x, 

J <r +x J 
a a 

x xx 

J G(py')' = Gpy'| - J G'.py' 

a a a 

_ x 

= p(x)y'(x) G(x, x ) - p(a)y'(a) G(x, a) - J (pG').y'; 

a 

x xx 

J (PG')y' = pG*.y| - J (pG')'y 

a a a 

_ x 

= P( x ) G'(x, x ) y(x) - p(a)G'(a)y(a) - J qGy, 

a 

from 1. Hence, 

x x 

- J Gf = J G[(py')' - qy] = -p(x)y(x) G'(x, x“) + p(x)y'(x) G(x, x~) 
a a 

+ P(a)[y(a)G'(x, a) - y'(a)G(x, a)] 

= -p(x)y(x)G'(x, x") + p(x)y'(x)G(x, x“) 
since the terms in the square bracket give 0 by 3, and the boundary 
conditions on y. 


b b 

Similarly, thinking of J = ^ J with ?/ > x, 

b b 

- J Gf = J G[(py')' - qy] = p(x)y(x)G'(x, x + ) - p(x)y'(x)G(x, x + ) 
x x 

and so 


b 

- J Gf = y(x)p(x)[G'(x,x + ) - G'(x,x~)] - p(x)y'(x)[G(x,x + ) - G(x,x")] 

a 

= y(x) 

since G(x, x ) =G(x, x ), by 2; and G'(x, x + ) - G'(x, x~) = - —by 4. 

P(x) 
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Note that the Green's function is symmetric. 


Problems for §7.8 

P.25, 26 are Exercises 7.8, #1, 2. 

Solutions to Self-Assessment Questions 

1. If y(x) = x™ is a solution then 

2 ,. m-2 m-1 m 

x .m(m - l).x + ax.mx + bx =0, 

1. e. 

m 

(m(m- 1) + am + b)x = 0 

and so m(m - 1) + am + b = 0. Conversely if m satisfies this quadratic 
equation then y(x) = x™ is a solution. 

2. Consider the argument at the top of p.283 and note that the Helmholtz 
equation differs from that of Laplace by the term Au. One then sees that 
the differential equation for f(p ) is that given. 


If x = A 1//2 p, y = p ±//2 f and y', y",... denote derivatives with respect to x; 
f', f",... denote derivatives with respect to p, then 


dp 


d 1/2 dp 
(p f) • 


dx 


1 - 1/2 1/2 - 1/2 
= (-p f + P t ) A , 


dy' dp 
dp * dx 

= (- y>~ 3/2 £ + o" 1 / 2 f' + t> 1/Z £“ )A _1 . 


Hence 


x 2 y" + xy' + x 2 - ((n + |) 2 )y = x p~ jp ~f + p 


2. 1 -3/2 -1/2., 1/2.-1 

- -o t + o t'+p r ) • 


).A 

+ a 1/2 p (|p" 1/2 f + p 1/2 f' )A -1/2 + (Ap 2 - (n + |) 2 )p 1/2 f 
= p 1//2 [p 2 f" + 2p f 1 + (Ap 2 - n(n + l))f] 

and so result follows. 

[Note: When proving A = B always consider whether it is simpler to start 
with A or to start with B.] 

3. If v = (z 2 - l) n , then 


v' = 2nz(z 2 - 1) n ^ 


so that 


(1 - z )v' + 2nzv = 0. 


Hence (1 - z )v" + 2(n - l)zv' + 2nv = 0. 
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Now differentiating n times and using Leibniz' rule, 

[(1 - z 2 )v (n+2) - 2nzv (n+1) - n(n - l)v (n) ] 

+ 2(n - l)[zv (n+1) + nv (n) ] + 2nv (n) = 0 

and therefore 

/(n+2) _ (n+1) , , 4 (n) 

(1 - z )v - 2zv + n(n + l)v K ' - 0, 

i. e. 

(1 - z 2 ) - 2z ^-[v* n *] + n ( n + !) v (n) = 0. 

dz 

Therefore and so —— = P satisfies the Legendre equation. 

„n . n 

2 n! 

°° k 

4. Assume that there is a solution of the form x(t) = 1 + Y a^t . Then 


or 


i .tel ,2 n k 

ka k fc = + l a k t (a o = 


((k + 1) a k+1 - a k )t k = t 2 . 


Hence 


(k + l)a k+1 - a k = 0 (k * 2) 


= 1 (k = 2), 


and so 


a = a = 1, 2a = a , 3a = a + 1, 4a = a ,... 

J. U Z x <J Z J 

Therefore, 

1 _ 3 ^ _ 3 

a 0 ~ a i lf a 2 2' a 3 J7J f a 4 " 2.3.4'* * * 

and consequently 

2 2 4 

X(t) = 1 + t + ~ + 3(|y + tj- + ...) 

= l + t+ ^+3[l+ I i + | T +| T + ... - (l + t+^)] 
= Be 1 - 2(1 + t + ^). 

Check: 


x(t) = 3e - 2(1 + t), 

t 2 + x(t) = t 2 + 3e fc - [2(1 + t) + t 2 ] = 3e fc - 2(1 + t). 
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5. In the form of Theorem 7.2.3 the equation is 

2 2 

w' = f(z, w) = z + w (t = z, x(t) = w(z)). 

2 2 

Since z + w is analytic in z, w for all z, w the result follows from 
Theorem 7.2.3. 


If one were to argue directly as in the 'Second Proof' on p.290 then we are 

2 2 

led to consider p , p^ so that the coefficients of z , w in the series for 
Mp i^> 2 

- -—-— ( | z | < | w I < p „) are bounded by 1. Now 

( p 1 “ z ) \P 2 ~ w ) 1 2 


MPi „ 2 


M 


(P x - z)(p 2 - z) (1 - z/ P± ) (1 - w/p 2 ) 


= M (1 + — + + . ..)(1 + — + + ...) 

P , 2 p 2 

1 Pi 2 p x 


z w z w zw . 

= M(1 + — + — + — T T + - + • • . ) 

P, P n 2 2 P P 0 

1 2 p. p _ 12 


M M 

and given any p , p > 0 we can choose M large enough so that —> 1. 

J- Z z z 

P 1 p 2 


Hence it follows that the solution of the problem is analytic for all t, 
i.e. it is entire in t. 


Solutions to Problems 

1. Differentiate Legendre's equation 

„ ,2 

2 d y dy 

(1 - x ) —=r - 2x + n(n + l)y = 0 
'2 dx 1 

d x 


m times with respect to x and put v 


d y 


dx 


and we obtain 


(1 - x 2 ) —^ - 2x(m + 1) ^ + (n - m)(n + m + l)v = 0. 
dx 

m / 0 

Now set w = (1 - x) v and express the preceding equation in terms of w to 
obtain 

2 2 
2 v d w „ dw r , , „ m 

(1 - x ) —- - 2x — + [n(n + 1)- -}w = 0. 


dx 


1 - x 
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2. Put x(t) = y (t) and consider first the equation 
x(t) = -f(t) (t > 0) 

with initial condition x(0) = b. From Theorem 7.2.1 with u(t, x) = —f(t) 
there is a unique solution 
t 

x ( t) = b - J f(r)dr. 

0 

Next consider the equation 
t 

y(t) = b - J f( t) dr (t > 0) 

0 


with initial condition y(0) = a. From Theorem 7.2.1 with 


u(t, y) = b - J f(r )dr there is a unique solution 
0 


y(t) = a + j (b - J f(u)du)dr 
0 0 

t T 

= a + bt - J (J f(u)du)dr. 

0 0 

Interchanging the order of integration in the repeated integral we find that 

t T t 

J (J f(u)du)dr = J (t - u)f(u)du 
0 0 0 

and hence result follows. 

3. Let w(z) = w(—) = W(? ). Then 

f 

j *, ^ 2 

w'(z) = W'(f ).^ = W'(? ).(- 2 } = - ? W(? )? 

z 

w" (z) = - ( 2? W' (? ) + ? 2 W" (<r )) ^ 

= ? 3 ( 2 W' (? ) + ?W"(f )). 

Put p(-) = P(?)f q(-) = 0(0/ r(-) = R(? ) and the equation becomes 

? f i 

? 3 (2W •(? ) + ?W"(f)) + P(? ) . (- ? 2 W' (f ) ) + Q(? ) W(f ) = R(? ), 


l. e. 


w»(? ) + ij - ^Li) w (n + 5111 w(? ) = 5111. 


? 




From the definition of ordinary point the preceding equation will have ? = 0 
as an ordinary point if and only if each of 
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2 P(?) Q(? ) R(?) 

? 2 ' 4 ' 4 

f ? ? 

are analytic at ? =0, i.e. ? = 
these functions. Equivalently, 


0 is a removable singularity of each of 
if and only if each of 


n 2 , v 
2z - z p(z) 

is analytic at z 


4 4 

z q(z), z r(z) 

= oo. This means that for large z, 


P(z) 


2 

- + 
z 


-2 


q (z) = 


-4 


-5 



z z 


where the power series in - converge for all large z. 


The equation for W(? ) has a solution for small ? , 
y>u ) = o Q + c 1? + c 2 t 2 + ... 

assuming that ? = 0 is an ordinary point. This is equivalent to saying that 

if z = oo is an ordinary point of the original equation then this equation 
has a solution 


C 1 C 2 

W( Z ) = c 0 + — + — + ... 

z 


which is valid for large z. 

4. The equation has solutions analytic in {|z| < l} from the earlier 


discussion in the section. Let w(z) = Y a^z and then 


(1 - z 2 ) ^ k(k - l)a k z k 2 - 2z ^ ka j < zk 1 + n (n + 1) £ a^z k = 0. 

The coefficient of z from the individual contributions on the left is 
equated to 0 gives us 

(k + 2)(k + 1)a _ - k(k - l)a, - 2ka + n (n + l)a =0 (k > 0) 
k + 2 k k k 

and so 

(k + 2) (k + Da k+2 = -[n(n + 1) - k(k + l)]a k 

= -(n - k)(n + k + l)a , 

K 
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i. e. 


k + 2 


(n - k)(n + k + 1) 

(k + 2)(k + 1) a k 


(k > 0). 


Two independent solutions are obtained by taking a^, a^ to be 0, 1 and 1, 0 
respectively. 


If n is a positive integer then n - k = 0 for k = n and then a a ., ... 

^ 3 n+2 n+4 

are all 0 and so we get a polynomial. 

Remark: If one assumes a power series solution then one sees what the form 

of this must be and one can check directly that there are two 
independent solutions analytic in {|z| < l}. 

5. On p.298 we see that 

z 

W(z) = W Q exp [- J p(?)d?]. 

z o 

If we continue W^(z), W 2 (z) round C to obtain u^(z), ^ 2 (z) then the effect 


on W(z) is found by continuing J p(? )d? round C. Denoting by 

z o 

Z Z 

[J p(? )df ] the value before continuation and [J p(? )d? ] the value 

z o 1 z o 2 

afterwards we obtain 

z z 

[J p(<r )d<r ] - [J p(Od? ] = J p(? )d? = 2n- is. 

z 0 2 z n 1C 

0 0 

Hence, in the notation of the question, 

z 


W(w x r « 2 ) = W Q exp(-[J p(? )d? ] ) 


0 


= W Q exp(-[ J p(f )d? ] - 27riS) 

2 o 1 


= W n exp(-[J P(a)d?] ). exp( -2n iS) 

1 


= W(w , w 2 )e 


0 

-2 ir iS 


6. If we set w(z) = w(-) = W(? ) as in the solution to P.3 then from that 


solution we see that 
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)W' (f ) + 


QU ) 


0 , 


W"(? ) + <- - 

f t 2 




where P(? ) = p(~), Q(f ) = q(-). The original equation will have a regular 
singular point at z = oo if the above equation has one at ? =0. Multiplying 
by f^ this will be the case if 2f - P(? ) has a zero at ? = 0 and ^ is 


analytic at f =0. Hence we require that for ? near ? = 0, 

2 2 2 
P(f ) = a j/ + a 2 ? + ...» Q(0 = b 2 ? + b 3 ^ + ... 


Therefore the original equation will have a regular singular point at z = oo 
if for large z, 



7. Consider Theorem 7.4.1 and the preceding discussion. We shall get 
formal solutions since the procedure stated will amount to solving 
zu" + f(z)u' + g(z)u = 0 

as in Theorem 7.4.1, where f, g depend on m = m^, m^, the roots of 


m(m - 1) + p(0)m + q(0) = 0. 

m i m 2 

If m^ - m 2 is not an integer, then z '<£( z), z «/> ( z ) , where <j>, are analytic 

at z = 0 and neither is identically 0, cannot have the same behaviour at 

M N 

z = 0. Note that if <p[z) = Az + ..., tp (z) = Bz + then M + m^ ^ N + 


8. Consider the equation 

w" + p(z)w' + q(z) = 0. 

This equation has z^ e C as an ordinary point if p(z), q(z) are analytic at 
z . Hence if every point z^ in the complex plane is an ordinary point then 
p(z), q(z) are analytic for all complex z, i.e. they are entire. 


From the solution to P.6 if z = oo is a regular singular point, then 
p(z) -* 0, q(z) 0 as z ^ oo. Hence p(z), q(z) are bounded and so constant 
by Liouville's theorem. Clearly the constant value in each case is 0 and so 
the equation reduces to w" = 0. 

[Remark: We use that if f(z) is entire and f(z) ->0 (z -» oo) then f(z) = 0. 
This is proved as follows. 


43 



To say f(z) ->0 (z -> o©) means that given <r > 0 there is an M such that 
|f(Z)| < e (|Z| > M). 

Since |f(z)| < s (|z| = M) it follows that |f(z)| < c (|z| < M) by the 
Maximum Modulus Principle (Theorem 3.7.1/ p.119). Hence f(z) is bounded and 
so it is constant by Liouville's theorem (Theorem 3.6.4, p.115). Clearly 
f(z) = 0.] 

9. When all points are ordinary apart from z = 0, oo, then 
w" + p(z)w' + q(z)w = 0, 

where p(z), q(z) are analytic for all complex z other than z = 0. If z = 0 

2 2 

is a regular singular point, then z p(z) has a zero at z = 0 and z q(z) is 
analytic at z = 0; and if z = oo is also a regular singular point, then zp(z) 
2 

and z q(z) are bounded near z = 0, by solution to P.6. Therefore zp(z), 

2 

z q(z) are entire and bounded and hence constant by Liouville's theorem. 

A B 

This means that p(z) = -, q(z) = —— for some constants A, B and then the 

z 2 

z 

equation is 
2 

z w" + Am + B = 0, 


If z = -, w(z) = W(? ) then the equation comes 


? W" + (2 - A)?W' + BW = 0, 

from solution to P.3. Hence the indicial equation at z = oo, f = 0 is 
n(n - 1) + (2 - A)n + B = 0. 

Therefore 

m, + in = 1 - A, n, + n„ = A - 1, 

12 12 


m i m 2 = B ' 


n i n 2 = B - 


and so 


m l + m 2 + n i + n 2 = 0, ryn., = n^. 


10. If all points are ordinary apart from 0, 1, oo which are regular singular 
then 

w" + p(z)w’ + q(z) = 0 

where p, q are analytic in the plane apart from z = 0, 1. When z = 0, 1 are 
regular singular then p(z) has at most simple poles at these points and q(z) 
has at most double poles at these points. Therefore 
, v A B 

P(z) = " + --y + <A(z) f 

z z - 1 


44 



C D E F 

g<z> =; + 7 + + 77777 


+ ^(z) r 


where <fi{z), ip{z) are entire. If z = oo is a regular singular point then 
9 

zp(z), z q(z) is bounded near z = oo and so are z<j>(z) r zip{z). Consequently 
<p(z), ip{z) are bounded and 4> {z) , ^(z) -> 0 (z ■> oo) so that <p{z) = 0, 

2 c E 

xbiz) = 0. Also z (- + -) is bounded near z = oo so that E = -C. Hence 

z z - 1 


the equation is 


A B . r D E 

W" + (- + -- W' + {—+ - 

z z - 1 2 

z ( z - 1) 


2 + C( 5 - F^T )} w ■ 


and this is of the form given in the question. 


Consider then the equation of the question. The discussion on p.303 and the 
form of the equation at z = 0, 1, « shows that 
m(m - 1) + a Q m + b Q = 0 has roots a 2 ? 

m(m - 1) + a ]L m + b = 0 has roots p p 2 ? 

and, using the equation transformed by z = -, 

m(m - 1) + (2 - a - a ± )m + (b Q + b 1 - c Q ) = 0 has roots y 7 2 * 


Hence 


“l + “2 = 1 ” V "l + P 2 = 1 ” a i' 
“l“2 = b 0' ^1^2 = V 

7 1 + 7 2 = a 0 + ®1 ' 11 7 l' , 2 = b 0 + b l 


and so 

c Q = b Q ’+ b x - 7 X 7 2 = + ^2 ~ 1 ^ 2 ' 

a l + a 2 + + 7 1 + 7 2 = 1 ~ a 0 + 1 ~ a i + a 0 + a i ~ 1 = 

P.ll. Take w(z) = £ r(a k) 77 in the left hand side of the 


equation and we get 


z(l-z) j> 
k = 2 


r{a+k)r(p+l(L) z 


k-2 


r ( 7 +k) 


(k-2)! 


OO 

[ 7 -(a+/? + l)z]££ 


r(a+k)r(/?+k) 
r( 7+k) 


k-1 

z 

(k-1)! 


— ap 


r(a+k)r(/?+k) 

r(7+k) 
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and the coefficient of 


k . 
z is 


A k = 


r(a+k+i)r{/3+k+i ) i 


r{a+k)r{/3+k) 


r{a+k+i)r{/)+k+i) l 


r( 7 +k+l) (k-l)! 

,T(a+k )T(/?+k ) 1 


- («+/? + !) 


r( T +k) (k-2)! 

r(a+k)r(f}+k) l 


r( 7 +k+i) 


r( 7 +k) (k-l)! 


a/3 


r( 7 +k) * k!* 


We set -y = 0 whenever n is a negative integer. Hence 
Ma+k )r(/? + 

r ( 7 +k) .k! 1 7 +k 


A. = r [« +k ) r (P +k ) [ (a+k)(/?+k) k _ k(k _ 1} + (*+k)( 0 +k) 


■+k 


(a+/3+l) .k - a/3 ] 


r(a+k)r( 0 +k) r («+k)( 0 +k) fl 2 

= — n^k) — [ ~— —~ (k+7) - {k + (a+/?)k + 


r{a+k)r( 0 +k) 


r( 7 +k) 


[(a + k)(/?+k) - (k+a)(k+/?)] 


= 0 


k! 


Hence F(a, 0, 7 , z) satisfies the equation. 


One can verify in the same tedious way that the other function is also a 
solution. 


The functions have factors r( 7 + k) and r(2 - 7 + k) and the condition k not 
an integer ensures these are non-zero, 

12. In the section it is proved that J (z), J (z) are linearly dependent 

v -v 

if and only if v is an integer. In this case W(z) = 0 and - ^ - in un = 0 

z 

and hence the result is true in this case. 


Assuming v is not an integer, first of all we argue as in the proof of 
Theorem 7.3.2 to obtain 


W( z) = 5 

for some constant A. Consider 


J_^(z) on p.312 (remember the 


n = 1) we see that 


j (z) . .. ( . v 2z > 

; r(l + i /) 


+ ..., 


ing the 
sum for 


J'(z) = 


initial terms of the series for J (z), 

v 

J (z) starts at n = 0 and not at 

V 


V (l/2z ) v 1 
2 r( 1 + v) + ' " ' 


and so 



W(z) = 


v/2 .(l/2z) 1 _ v/2 (l/2z) 1 + 

r(i + i/)r(i - i/) l r(i - i/)r(i + i/) J 

2v 1 

r(i + i/ )r(i - v) *z + “ *' 


and hence 


W(z) 


2v 1 

r( l + v )r(i - i/) 'z 


On p.197 it is proved that 


2 

r{v )T(1 - £/)* 


r(z)T(l - z) = n cosec n z 


and 


taking z 
W(z) = 


= v we 
2 sin 

n- Z 


obtain 

UTT 


n 

sin n z 


13. Putting t/ = ^ in the sum for J^(z) we get 

00 , , . n ., . 2n+l/2 

t t 7 \ = V (-D d/2z) 

1/2 1 ; A n!T(1/2 + n + 1) 
n=0 

and so 


,z.l/2 _ . (-l) n (l/2z) 2ntl 

( 2' 1 / 2 ' ' i. n!T(n + 3/2) - 


n=0 


Now 


3 4 1 , 1, 3 1 „.l. 

P(n + ^ = ( n + 2 * n 2 * * * 2 * 2 ^ 2 


(2n + 1)(2n - 1).. .3.1 
»n+l 


■Jn f 


putting z = | in the result on p.197, and 


n! = n.(n - 1)....1 

2n(2n - 2) . . .2.1 


and hence 


3 X (2n + 1)! , 

nlr(n + 2 = n+1 ' / " 


Using this result in the sum above we obtain 

1 


n 2n+l 

(-1) Z sin z 


,Z 1/2 -L y 

2 J i/ 2 (z) = 77 J 0 (2n + ryr = -v? 

and therefore 

2 1/2 

J, /0 (z) = ( —) sin z. 

1/2 n Z 

The result for J^^fz) is P roved similarly. 
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N | l-> 



I 

Since ^ is not an integer the solution of the given differential equation is 


/ c sin z + c cos z 

Vl/2< Z > + C 2 J -1/2 (Z) = -7i- 

A sin z + B cos z 
\/z 


14. Taking v - | in the identity 

J ,(z> + J +1 (z) = “ J (z) 

y — 1 UTJL 7. V 

on p.315 we get 

J 3/2 <Z> = z J 1/2 (Z) " J -1/2 (Z) 

= ,^)V 2 ["A"-* - COS z], 
7T Z Z 

J-3/2( Z ) iS dealt with similarly. 

15. We have that 


J <z) = l ( ~ 1)n(Z / 2 - ^ - = 1 - (?) 2 + 


n=0 (n!) 


and so 


V z) = ; J 0 (z) I —— df ' 


? [j Q (f )] 


as on p.313, has a term - J n (z).log z and so the result follows 

n U 

16. Note that the equation should read 


(xy' )' y+Axy=0 


or 

2 2 2 
x y" + xy' + (Ax - n )y = 0. 

Put Vax = y(x) = y(^-) = Y(€ ) and then 

y • (x) = Va Y'((), y" (x) = AY" U ) . 
Consequently 

2 e £ 2 2 

— .AY" + -r-.y/\Y' + (A — - n )Y = 0, 
A VA A 

i. e. 

2 2 2 
£ Y" + (, Y 1 + U - n )Y = 0. 
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This is Bessel's equation of order n and from the discussion on p.313 the 
only solution bounded near the origin is a multiple of 
J n U) = J (>/Ax) when n is an integer. 


For the final part of the problem put 

y.(x) = J (x/A . x), y. (x) = J (Jx . x) 
i n 1 3 n j 

and consider 


(xy!)' - — y. + A. xy. =0 

1 X I 11 


(xy ) '-y . + A . xy . = 0 . 

3 x 3 3 3 

Multiply the first equation by y. and the second by y, and subtract to 

3 l 

obtain 

(X yp'Yj - (X yll'y. = (A ^ - A.)xy.y . 

Now, integrating by parts, 

b b b 

J (X y! ) 'y. = x y! y. | - J x y! y'. 

0 3 3 0 0 1 3 

b 

= "J x ^i 

0 3 

since the integrated term is 0 at x = 0 because of the factor x and is 0 at 

x = b since y .(b) = J (Va . b) = 0. We have the same result with i, j 
3 n ;j J 

interchanged and so 

b 


(A . — A . ) P x y. y , = 0 
3 1 l 1 


0 

and if A ^ ^ A the desired result follows. 

17. Put t = i<9 and we get 

1 r z s i- n h t - nt 1 / z sinh i0 - im9 

571 .J e dt = 27 f e d * • 

— llT —n 


Since 


ie _ -i# 

sinh i(9 = -—-- = i sin e 


we obtain 


1 1?r 

iTT I 

— ITT 


z sinh t - nt ... 
e dt 


1 

2jt" / 


i ( z sim9 - ne ) 


d# 
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= j (Z), 

n 

as shown on p.316. 

18. Consider 

(1 - x 2 )P" - 2x P' + n(n + 1)P = 0, 

n n n 

(1 - x 2 )P" - 2x P' + m(m + 1)P = 0. 

m m m 

Multiply the first by P and the second by P and subtract to obtain 

m n 


(n(n + 1) - m(m + 1) f P P = f (1 - x)(P" P - P P") 

J n m J ' m n m n 7 


-1 


-1 


+ 2 r x(p' p - p' p ) 

n m m n 


Now 


1 2 2 1 1 2 
f (1 - x )P" P = (1 - x )P' P I - f (- 2x P + (1 - X )P')P ' 
J _ 1 ran m n _ J ' n n 7 m 

1 1 2 
= 2 r x p p* - r (i - x )p- p* f 

J ^ n m J nm 

with a similar result with m, n interchanged. Hence 

1 

(n(n + 1) - m(m + l))f P P = 0 

J ^ n m 

and when m ^ n the result follows. 

19. Consider Rodrigues' formula 

1 d r 2 n, 

P (x) = —-—[ (x - 1 ]. 

n „n . , n 

2 n! dx 

Then 


, ry 1 .H 

J- . 2* /• d 


_ n 

r n v-, , i .2 r c z n. d 2 . v n._ 

P (x) dx = (——) | —-[ (x - 1) ].——[ (x - 1) ] dx 

J i n J , , n _ n 

-1 2 n! -1 dx dx 

If we integrate by parts n times the integrated term will be 0 at 


x = ±1 at each stage because of a factor x - 1 and each integration will 

introduce a factor -1. Hence 

1 , , ,n 1 ,2n 

r 2 _ (-1) r 2 n d 2 _ .n. 

j P M ( x )dx = —-- - | (x - 1) .——[ (x - 1 ]dx 

-1 (2 n!) -1 dx 

and 


,2n _ 2n 

2 ,n, d r 2n , ^ 

-[(x - 1) ] = —— [x + lower powers of x] = (2n)! 


dx 


2n 


dx 


2n 
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and so 


. 2 _ ,n , _.n,, 2 . n 

, using (x - 1) = (-1) (1 - x ) , 

f P 2 (x)dx = - - - ( - 2n) 1 f (1 - x 2 ) n dx 
-1 n ( 2 n n!) 2 -1 

2 . (2n) ! p 1 2 n 

= —--~ I (1 - x ) dx, 

( 2 V ) 2 0 

since the integrand is an even function. Put x = cos 6 and then 

1 9 n / 2 o n4 .i 

J (1 - x ) dx = J sin 6 dtf. 

0 0 


We have 


r 7r/2 . 2k+l / /2 . 2k-l r /2 . 2k-l 

sm 6 d(9 = J sm 6 66 - \ - 

0 0 0 


sin 


6 cos 6 66 


and 


n/2 , . 2k ir/2 7T /2 . 2k + l 

r • 2k_ l„ „ „ sm 6 1 r sm 6 

j sm 6 .COS 6 .COS 6 d0 = -—-. COS0 I + 

0 2k 0 


• J 


2k 

n/2 . 2k+l 


0 


2k 


sm 


Therefore 

n /2 


• 2k+1 „ ^ 

sin 6 66 - 


2k f 

2k + 1 J 


0 


n /2 


2k 


66 . 


. 2k-l „ 
sm 6 66 . 


0 0 
Repeated applications of this result gives 


r ?r/2 . 2n+l 2n 2n - 2 

i sin 6 66 = 

0 


2n + 1 * 2n - 1 .3 * 


H 


n/2 


sin 6 66 


[2n.(2n - 2) ... 2] 

(2n + 1).2n(2n - 1) ... 3.2 

( 2 n n !) 2 
(2n + 1)!* 


From the above 
1 


r P 2 (x) dx = 2 -< 2n > ! ( 2rln! > 2 

-1 ?n ( 2 n n!) 2 (2n + 1)1 


2n + 1 

20. We use Rodrigues' formula 


, ,n 9 

P (x) = —-- —[(x - 1) ]. 

n n . , n 

2 n! dx 


66 
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The polynomial (x 2 - l) n has n zeros at x = -1 and n zeros at x = 1. By 
(3 2 n 

Rolle's theorem — [(x - 1) ] has n - 1 zeros at x = -1, n - 1 zeros at 

dx 

x = 1 and a zero x^ 1 * satisfying -1 < x^ < 1. Again applying Rolle's 


d ^ 2 n 

theorem, —-[(x - 1) ] has n - 2 zeros at x = -1, 

dx 2 

n - 2 zeros at x = 1 and a zero x| 2 ^ satisfying -1 < x| 2 ^ < x^^ and a zero 

*2 ^ satisfying x| ^ < x^^ < 1. We continue this argument for n steps and 

note that at each stage the zeros in the open interval (-1, 1) are distinct. 
At the nth stage we see that 
-.n 

d _ , 2 „ 4 n, t 

-[(x - 1) ] has n distinct, i.e. simple, zeros in (-1, 1). 


dx 11 

21. If we start with the associated Legendre equation instead of the 
Legendre equation as we did in the solution to p.18, then the argument 
proceeds in precisely the same way as in the earlier solution. 

22. We argue as in the solution to P.19. On p.323, 

T-% m , . , . x n /n 2 m/2 d m „ , v 

P (X) = (-1) (1 - X ) -- P X) 

n , m n 

dx 

_ _m+n 

= (-1 ) n (l - x 2 ) m/2 .-J- -1—[( X 2 - l) n ] . 

n m+n 
2 n! dx 

2 n 

If m > n, then m+n exceeds 2n, the degre® of (x - 1) , and so the above 
derivative is 0. It is usual, cf. the /"-function, to put 
1 


(—k ) 1 =oo, 


0 for k = 1, 2, ... . Hence the result of the problem is 


(-k)l 

true for m > n and so we assume that m < n. 


Consider 

1 

J [P (x)] dx 
-1 

, _ 1 _ .m+n _ _,m+n 

- ( ^r~> I C(1 - x > m ' 1 > n)) - J ^ [(x - 1)n]dx - 

2 n! -1 dx dx 

We integrate by parts m+n times with the differentiations applying to the 
factor that arises form the function in brace brackets. This function is a 
polynomial of degree 


52 



2n - (m + n) + 2m = m + n 
with highest coefficient 


m 


(-1) 2n(2n - 1) ... (2n - m - n + 1) 

and so after m + n differentiations we arrive at the constant 


m 


(m + n)!.(—1) 2n(2n - 1) ... (n - m + 1) 

- (-l) m .(2n,. < n + m)! 


(n - m)! 

At each integration by parts the integrated term is 0 at x = 1, 

2 

x = -1 because of a factor x - 1 and a factor (-1) is introduced. Hence 


/[p”( x )] 2 d x . <^)Vl> m (2n,i J <** - 1) 

-1 2 n! -1 


(n + m)! m+n 


dx 


(2n)! (n + m)I - 2 n 

, 7~ -— TT J (1 - X ) dx 

(2 n n!) 2 (n - m)! -1 
2 (n + m)! 


2n + 1 (n - m)1 


using the value of J (1 ~ x 2 ) n dx found in the solution to P.19. 

-1 

23. The left hand side of the given equation is of the general form on p.325 
with 

— 3 x — 3 x 

p(x)=e , q(x) = 0, p(x) = e 

Considering the transformation at the top of p.326, where 

X 1 /n /2x 

t = f (-) x dx f we can take t = x and u(t) = e y(x). Then 

J p 

u(t) = - | e 3//2x y(x) + e 3//2x y'(x), 

u(t) = 9 - e" 3/2 y(x) - 3 _3/2x y'(x) + e“ 3/2x y"(x) 

= | u (t) + e 3//2x (e 3x y' (x) ) ' 

= | u (t) + e 3//2x [f(x) - Ae 3//2x u(t)] 

and so 

9 . 3/21 _ /1- , 

u + (A - -)u = e f(t). 

9 3/2t 

Therefore q(t) = - q(t) = e J f(t). 

24. Here we are considering a = 0, b = n and a. = /) = 0 in the context of 
(a) f (b) on p.326. The homogeneous equation y" + Ay = 0 has the general 
solution 


53 



-AA sin Ax + AB cos Ax 


y(x) = A cos Ax + B sin Ax, y'(x) = 
and hence 


,, . sin Ax , . 1 . r ^ 

<Mx, A) = —--, ^(x, a) = - sin iA (tt - x)j. 

A A 

<P f A are obtained by choosing A, B so that y(0) = 0, y'(0) = -1 and 
y(/r ) = 0, y' (tt ) = -1 respectively as on p.327. In this case 

W {(p , xp) = <pxp ' — <p ' xp 


sin Ax % . . . sin A (* - x), 

---.-cos A (?r - X) - [- COS AX.----] 


= -[sin Ax cos A (n - x) + cos Ax. sin A (tt - x)] 

A 


sin A tt 
_ 


W(A). 


Hence the eigenvalues are ±1, ±2, ... . 


If A = k, a non-zero integer, then 

k-1 

, . iv sin kx , , 1 (-1) 

0(x, k) = - -^—, ^(x, k) = - sin k(^ - x) = --- sm kx 

k—1 

and so, in the notation of p.332, A = (-1) . Also 

K 


w' (k ) 
# k (x) 


k—1 

-cos kn = (-1) and consequently, as on p.333, 
. A k ,1/2 , 

'^TET 1 Hx ' k) 

sin kx 


From this the form of the expansion result in this particular case follows 
from Theorem 7.7.4. 

25. One can check that A = 0 is not an eigenvalue of the Sturm-Liouville 
problem and so the Green's function is determined by properties 1-4. on 
p.340. In this case p(x) = 1, q(x) = 0 and so 1. gives, for each x in 
(0, 1), 

At + B (0 < t < x) 

G(t) = G(x, t) = 

A 2 t + B 2 (0 < t < 1), 


where A^, ..., B 2 are independent of t. Property 2. shows that 


Ax + B = Ax + B_ ? 
112 2 

and property 3. shows that 
B 1 = °- A 2 + B 2 = ° 
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since in this case a = 0, b = 1 and a = p = 0. Finally property 4. shows 
that 

A 2 ' \ ■ - 1 ' 

Solving for A , . .., Bwe obtain 

(1 - x)t (0 < t < x) 

G(x, t) = 

x(l-t) (x < t < 1). 


For this problem p(x) = 1 and so, 
K(x, t) = x/p(x)p(t) G(x, t) 
u(x) = Vp(x).y(x) = 


as on p.342, 
= G(x, t), 
y(x)r 


1 

F (x) = J [Jp(x)p(t) G(x, t) .f {t)/y/p (t) ]dt 
0 

1 

= J G(x, t ) f ( t )dt, 

0 


and hence 

1 

U(x) = F(x) + A J K(x, t)u(t)dt 
0 


is the associated Fredholm integral equation. 


The kernel K is determined from the homogeneous equation and F is determined 
from f and the homogeneous equation. The solution y(x) of the differential 
equation with boundary conditions and the solution of the integral equation 
are related by u(x) = v/p(x) y(x) = y(x) so the solution of the problem in 
one form gives immediately the solution in the other form. 

26. One can argue directly as in the earlier problem or first obtain an 

f , £ 

equivalent problem as follows. Put x = e , (, = log x and then, with y(e ) = 
YU ), 

y(x) = Y'u >||, xy' (x) = S' ({ ); 

g|[xy'(x)] = *"<«>§! = e ’ 4 y"(f)• 

Hence the equation of the problem becomes 

e _ * Y"U ) + Ae^Y(0 = -f(eS (0 < ^ < 1) 

or 

Y"(€) + Ae 2 *YU) = -e^f(e^) (0 < i < 1), 
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with boundary conditions Y(0) = Y(l) = 0. The Green's function in this case 
is the same as in the previous problem, but now 

p = p{n = e 2 * and f is replaced by e*f(e*)- One can now find the Fredholm 
equation in this case for U(?)r as on p.342, and setting U(log x) = u(x) and 
changing the variable of integration from t on [0, 1] to r on [1, e] by the 

transformation r = e we find the Fredholm equation for u. 
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Chapter 8 


Fourier Transforms 


| READ §8.1 | 


This section is rather well presented and deals with a number of the basic 
results for Fourie,r series. 

Self-Assessment Questions 

SAQ 1, 2, 3 are Exercises 8.1, #1, 2, 3. 


| READ §8.2 | 


Theorem 8.2.1 The 'Cauchy principal value' means that we consider 

lim 1 r r 

R-*x> 2?T J J as in the fourth line of the proof. For the integral to 

—R —oo 


R. 


lim 


exist in the ordinary sense we would require J "“‘ l — f f where 

R /R-^oo 2n J J 

1 2 —«! —oo 


we 


consider R^, R 2 ^ 00 independently. 


dx 


For example, J -- exists in the latter sense since 


-oo 1 + x 


R„ R 

2 2 
r dx ^ -1 - -1 -1 

I -- = tan x I = tan R„ + tan R. 


-R 1 + x 


-R. 


n 


as R l f r 2 * 00 in an Y manner, i.e. given e > 0 there is an R (c) such that 


for all R x , R 2 > R Q ( £ ), 




dx 


— 7T < e . 


-R 1 + x 


On the other hand for R,, R > 0, 

12 
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xdx 


xdx 


xdx 


r xdx r 

J _ " 2 ~ J 


-R 1 1 + x 2 J 0 1 + x 2 J 0 1 + x 2 


1 1 + R ? 

= 5 lQ g(-2 J 

1 + R. 


xdx 


and so -~ does not exist m the usual sense* However, if R = R , 

-oo 1 + X 2 12 

.. xdx 00 vdx 

then j - - = 0 so that the Cauchy principal value of J -- is 0. 


-R 1 + x 


-OO 1 + X 


Self-Assessment Questions 

00 

SAQ 4. If J | f(t) | dt exists prove that G(x) = 1 


J 2n 


J f (r )e 1X7 dr is 


continuous on (-oo, oo) and G(x) 0 as x -* ±oo . 
SAQ 5. Exercises 8.2, #1. 


| READ §8.3 | 


Theorem 8.3.4 This result should be compared with Theorem 8.1.3. 

Problems for §8.1, §8.2, §8.3 
P.1, 2 are Exercises 8.1, #6, 7. 

P.3, 4 are Exercises 8.2, #2, 3 

P.5, 6, 7, 8, 9 are Exercises 8.3, #1, 2, 4, 5, 6. 

[Note: Exercises 8.3, #5 is badly stated. Part of P.8 is finding the 
correct statement of this exercise.] 

< ♦ 


| READ §8.4 | 


Theorem 8.4.6 At the start of the proof there is differentiation under the 
integral sign. When the limits of integration are oo and —oo justification of 
this can be quite tricky. Generally speaking if this operation leads to a 
convergent integral, then one expects this operation to be valid. 


58 



is 


7 t 

p.372 The argument used to show that e h(t) is piecewise continuous 
rather unsatisfactory. It is best to assume that those f, g to be 
considered will always lead to an h with sufficiently nice properties for 
the application that is being considered. 

Problems for §8.4 

P.10, 11, 12, 13, 14, 15 are Exercises 8.4, #1, 2, 3, 4, 5, 6. 


| READ §8.5 | 


Theorem 8.5.1 In this theorem we consider 

-l.„ r „ -1, F[f]. 

y = F [F[y] ] = F [ p ( - j z -y] 

in a somewhat formal manner under hypotheses that ensure that the integral 
for the inverse on the right exists. But instead of proving this integral 
is the inverse, the result is verified by assuming differentiation under the 
integral sign is permissible. A proper justification for such an argument 
can be rather tricky in a general context and it is perhaps better to check 
that everything is all right in each particular case. 

p. 379 The discussion on the one-sided Fourier transform is to be superceded 
by the Laplace transform in the next chapter and so one should not spend too 
much time on the details. 

Problems for §8.5 

P.16, 17, 18 are Exercise 8.5, #1, 2, 3. 


| READ §8.6 | 


The examples in this section illustrate the way in which the theory of 
Fourier transforms can be applied to partial differential equations. Note 
that the emphasis is on Fourier transforms and not on the 'easiest' or 
'best' way of solving a particular problem. For the wave equation, for 
example, d'Alembert's solution should always be kept in mind, i.e. 

2 2 

, . £ a u l a u 

u(x, t) = f(x + at) + f(x - at) is a solution of —- = —j —y* 

ax a a t 
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The solutions of the examples could be better presented in some cases and 
various implicit assumptions are made. It is better, I think, to regard the 
arguments given as routines for arriving at possible answers which require 
to be checked directly. However, in some cases such verification can be 
rather difficult. 


As an illustration of the points made above consider Example 8.6.1. To 
obtain 

2 1 d 2 U 

(iz) U(z, t) = — —^ 


a d t 


we consider 


r d U -1ZX, 

J —- e dx 


00 2 

1 r d U -izx n 

— f —- e dx. 

2 J 2 

a — oo a t 


—oo d x 

The integral on the left is integrated twice by parts to give 

00 2 , 00 00 

d —U -izx = £U e -lZX ! + iz r LH. e - lzx dx 

J 9 d v j 


-00 d X 


dx 


, OO . w 0 

— e 1ZX | + iz.ue 1ZX | + (iz) J ue 

ax J . 


-1ZX 


dx. 


The given condition on u(x, t) as |x| oo ensures that the second term on 
the right is 0 when z is real, but not the first term. For this we require 

d U 

an additional condition, viz. — ** 0 (l*| * 00 ) • 

a X 


To obtain the equation for U we must assume that 
0 
J 


2 2 oo 

d u -izx, d r -izx 

—_ e dx = —- J ue dx 

—oo d t dt -oo 


a 2 u 


d t 


i.e. that the operations of integration from -oo to oo and partial 
differentiation of the 2nd order can be interchanged. This would require 
quite severe a priori restrictions on the 'unknown' function u(x, t) if one 
wished to be rigorous. 


Example 8.6.1 illustrates some of the main features of the use of Fourier 
transforms in §8.6. One has two variables x, t and one takes the Fourier 
transform of u(x, t) with respect to one of them. After integration by 
parts and an interchange of differentiation and integration one obtains an 
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ordinary differential equation. In each integration by parts the integrated 
part must be 0 at both limits. in order to ensure this, the transformation 
variable must be suitably chosen and an appropriate form of the transform 
used. 


The ordinary differential equation is solved taking account of the initial 
and boundary conditions and the solution of the original problem is then 
found by applying the inverse Fourier transform. It may not be clear to 
begin with which of x and t should be chosen as the transformation variable. 
If, having chosen one of x and t, one runs into an impasse then one must 
start again with the other variable. 


It may happen, for example, that for t > 0 one can check that a 'solution' 
one has obtained is correct, but at t = 0 the solution has some kind of 
peculiar singularity. In some way this may reflect a peculiarity of the 
physical system. When one shakes a string or applies heat the physical 
phenomena at the 'interface' are probably extremely complex and perhaps a 
strange singularity of a solution is the mathematical counterpart of this 
complexity. 

Problems for s 8.6 

P.19, 20, 21 are Exercises 8.6, #1, 2, 3. 


§8.7 is not part of M828. 


Solutions to Self-Assessment Questions 


1. Define F(t) - f ( t) so that F on( -it , n) corresponds to f on (-L, L), 


i.e. F(?r) = f(L), F(-tt) = f (—L) . Applying Theorem 8.1.5 we obtai 


m 


F(r + ) + F (r ) A 0 


2 ■ \ (& k COS k r + B sin k r ), 

k=l 


where 


1 * 1 n 
= - J F(t) cos kt dt, B. = - f F(t) sin kt dt. 
it J K nJ 


Put u = — and then 

n 

L 

a 1 r *, . k* u 

\ ; J f(u) cos — 

-L 


_ 

du = - J f(u) cos du = a , 
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with a similar result showing that B, = b, . Finally put t = -r in the 

k k L 

series above to obtain the desired result. 

2. In the series of the preceding problem put 

. ikrr /Lt — ik;r /Lt 

k?i- e + e 

cos — t = - 2 -< 

. ikr/Lt -ikrr/Lt 

kn- ± e - e 

sin — t = - —- , 

L 21 

to obtain the series of this problem with 
_ a k b k _ 1 

C k 2 + 2i 2 a k lb k ' 

a i b i 

_ k k 1, .. 

C -k = ~2 “ JI = 2 U k + lb k } * 

3. We have the result of the first problem with 

1 L kn 

a = - f f(t) cos — t dt 

k L J L 

~L 

1 L ° 

= f (f + J ) f (t) COS t dt 

L 0 -L L 

and if t = -r , 

0 kr ° k* 

f ( t ) COS — — t dt = J f (—r ) COS ( —. -r ) (-dr ) 

-L L 


f f (r ) COS — n —.r dr 

o L 


Hence a^ = 0 and similarly 

1 L kn 1 L 0 kn 

f(t) sin —— t dt = - (J + J ) f(t) sin — t dt 
-L 0 -L 

1 L L k, 

= r (J + J ) f(t) sin — t dt. 

0 0 

4. Given s > 0 choose X = X(c) > 0 so that 
-X 00 

J | f ( t ) | d t < £ , J | f ( t) | dt < e. 

-oo X 

00 -X X oo 

Then, considering J = J + J + J , 

—oo —oo —X X 


|G(X + h) - G(x) | < -jj- + -7^-11 f (r ) [e l(x+h)7 _ e 1Xr ]dr 
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4e K max ■ -ihr . 

s 75T + 75T -x<r<x |e ' 11 ' 

00 

where K = J |f(t)|dt. Now e is continuous in t and is 1 when t = 0 and 

—00 

so there is a t^ = > 0 such that 

|e lt -l|<e (|t| < t Q ) 
and hence 

_ ^ t 

Ie r - l| < e (~X < r < X, Ih! < — = h^, say). 

Consequently, 

|G(x + h) - G(x)| < -jj- (4 + K) 
for |h| < h , i.e. G is uniformly continuous on (~oo, oo). 


From Theorem 8.1.1, 

X . X X 

— i Xr 

I f(r)e dr = J f (r ) COS Xr dr + i J f(r ) sin Xr dr 
-X -X -X 

■+ 0 (X -> oo) 

and so 

X 

| J f (r ) e 1Xr dr I < £ ( I X I > X ) , 

-X 

where x^ = ) > 0. Hence 

00 

| J f (r ) e 1Xr dr I < 3e ( |x| > X Q ) , 

—OO 

i.e. G(x) + 0 (x ±oo). 

Remark If a < b are finite and f is Riemann integrable over [a, b] then 
given e > 0 there is a function g continuous on [a, b] such that 
b b 

| J f(t) dt - J g(t) dt| < e . 
a a 

This result shows that Theorem 8.1.1 is true for any integrable function. 
5. The Fourier transform is 

00 

_ . lr “3r — lXr _ 

G (x) - ~^2 n J e * e 
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1 r -(a+ix)r , 

757 l e dr 


Jin * a+ix 


r~. e 


-(a+ix)r 


1 _1_ 

Jin’ a+ix * 


izt 


For t > 0 consider the integral of G(z)e round the contour C of Figure 
8.2.1 on p.363. As R -► «>, the integral over the semi-circle tends to 0 and 
therefore 

R 

-i— I 

R->oo 


lim 1 ‘ v ixt, 2n i 

757 J_ R G(x)e dx ' T57' R ' 


izt 


where R is the sum of the residues of — r— in the upper half-plane. This 

a+iz 

1 —at 

function has a simple role at z = ai and so R = 76 and 

00 , 

1 r r./ v lxt ~ at 

j G(x)e dx = e 


For t < 0 we consider the semi-circle from R to -R in the lower half plane 

R 

and for t = 0 we evaluate J G(x)dx directly and then let R <*>. 

-R 


Solutions to Problems 

1. If f'(t) has the Fourier series 

A oo 

—- + Y (A. cos kt + B. sin kt) 

2 . L v k k 

k=l 

then, if t is a point of continuity of f'(t), 


f ' ( V = ~T + I cos kt Q + B k sin kt Q ), 

k = l 

by Theorem 8.1.5. Now 

TT 

1 


A. 


- f f ' (t) cos kt dt 

n J 


1 *■ k n 

= -f(t) cos kt I + - f f(t) sin kt dt 

n TT J 


= k -V 
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where is the kth Fourier sine coefficient of f(t). The integrated term 
is 0 since f (n- ) = f(-rr) by the periodicity of f(t). Similarly 

B k ■ ~ k v 

where a^ is the kth Fourier cosine coefficient of f(t). Also 


A 


o = ;J f,(t > dt 


= - [f(.) - £(-»)) 

n 

= 0. 


Hence 


■to = | 


[k b cos kt - k a, sin kt] 


and this series is the Fourier series of f(t) differentiated term by term. 
2. If a periodic function is differentiable so is its derivative since 
differentiating f(x + 2n ) = f(x) gives f'(x + 2n ) = f'(x). Now 


a k 


= - r f(t) cos 

ir J 


kt dt 


1 sin kt 


f(t> i 


krr 


j f ' (t) sin kt dt 


■kT'J f ' (t) 


sin kt dt. 


Repeating this procedure n times we find that 

a k = ± -i- J f (n) (t) [?] dt, 

k n —n 

where [?] is either sin kt or cos kt. From Theorem 8.1.1 it follows that 

k n a ->-0 (k -► oo) and k n b is dealt with in the same way. 

K K 

3. Note that 


, -t 2 /2, 1 r ° 

F[e ] = 727 J 


-t /2 -ixt_. 
e .e dt 


-x 2 /2 


“72 n J* 


e"l/2(t + ix) dt# 


Z 

Suppose, for convenience, that x > 0 and consider Je ^ (?+ix) ^ round the 

pictured contour with ? = t, + ir/ . 
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>7 




_<- 

-R-ix 


R-ix 


The integrand is an entire function of f and so the integral is 0, The 
integrals over the segments [R, R - ix] and [-R - ix, -R] tend to 0 as R -> oo 
and so 


-1/2(t+ix) 


oo-ix 


at = J 


-oo- IX 


e -l/2(?+ix) 


2 

d? 


00 



e 





df 


—00 


- \J2ir . 

For the value of the final integral note that 

00 -t 

^ = r( 2> = J o dt 

oo 2 , 

= J e_{ /2 -« d < (t = « 2 /2) 

00 2 

= V2 J e _ * /2 d<; 

0 

and so 

00 2 00 2 . 

J e~* /2 d* = 2 J e" c /2 d$ = 2.V|. 

-oo 0 

-t 2 /2 

Hence result for F[e 7 ]. 

4. If f is real-valued then we can replace the integrand in Theorem 8.2.1 

by its real part and 

_ rc/ . ix(t-r). v „ . ix(t—r) 

Re[f(r)e ] = f(r ) Re[e ] 

= f(r ) COS[x(t - r )]. * 

Note that 
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f(r) COS[x(t - r)] = f(r)[cOS Xt.COS Xr + sin Xt.sin Xr ] 
and if f(r) is an even function then for x, t fixed f(r).sin xt.sin xr is an 
odd function. Hence the integral of this function with respect to r over 
(-oo, oo) is 0 and this is the second of the results asked for in the notation 
of Theorem 8.2.1. The last result is proved similiarly. 

5. Note that 

i \ m r | t | ^ 1 f* — | r | — 1 r Z _ 

G(z) = F[e 1 ] = -jj- J e 1 1 . e dr 


If z = x + iy, then 

Re [ - | t | - ir z ] = 

where 

0(r = 0) 

a - 1 (r > 0) 

-l(r < 0). 



—r — ir Z _ 
e dr 



r -ir Z 
e 


dr ]. 


- | r j + r y 

- | r | [1 - <ry] , 


The 

1 - 

For 


00 - Ir I - ir Z 

integral J e 7 dr will converge at 


ay > 0, i.e. provided |y| < 1 and so G(z) 
{z: z = x + iy, -oo < x < oo, -1 < y < l}. 
such z, 


J 


-r ( 1+iz ) 


dr 


1 + iz 
1 


1 -r(l+iz) 
t—. e 


1 + iz' 


r (1-iz) 


dr = 


1 - iz 


1 r(1 - iz) 
e 


both limits provided 
is analytic in the strip 


1 - iz * 

For the final part of the problem we must consider for -1 < y < 1, 

oo+ i 7 , j ^ oo+ i 7 


1 r itz o/ x J 1 r ' itz ,2 X 

727 J j _. e G(z) dz = 757 I . e -2 

—00+1 7 —00+1 7 1 + Z 


dz 


, oo+ i 7 itz 
1 r 6 


± J 

7T J 


dz 


-00+ i 7 1 + z 

From Corollary 8.3.1 we take 7 = 0 for convenience. 
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i tz 


For t > 0, consider J -— dx round the contour of Figure 8.2.1. For 

1 + z 

R > 1 the integrand has a simple pole at z = i with residue 


it z 


-t 


2 z 1 . 2 i* 

z = 1 

As R ■* oo, the integral along the semi-circle tends to 0 and so, letting 
R 00, 

00 itz -t 

c e , „ . e -t 

J —2 dZ = 2,rl '~2i = ' e ' 

-00 1 + Z 

and so the inverse in this case is 

1 „ _t -iti 

— .n e = e = e 


For t < 0 we proceed as above with the contour having a semi-circle in the 
lower half-plane. 

6 . In this case 

00 

j 11Z 

G(z) = F[u(t) cos wt] = -t^— f cos wt.e dt 

' /2 ' o 

, oo iw t —iw t 

1 r e + e ~* 1 1 z _ 

■ 757'^ -2-- e dt - 

As in Example 8.3.1, p.366, 

G < z > =2757 1 1 1 ’ 


w — Z w + Z 


1 2 iz 

2J2n'~2 2 

u — Z 

1 iz 

727* ”2 2 * 

u — Z 

To obtain this result we require 

i(w-z)t _ -i(w+z)t . v . 

e ■> 0, e ■> 0 (t ^ oo) 

and so we require Re[iz] > 0, Im z < 0. 


The inversion of the transform is essentially the same as in Example 8.3.1 
7. We are considering 

oo+ i 7 


, oo+ i 7 izt 

1 r e 


1 r , v izt , j. r 

757 J . G(z) e dz “ 57T I . 

—00+I7 — 00+I7 


dz 


in the two cases: 7 > 0 and 7 < 0. If 7 > 0, then 
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00+1.7 izt 00 ixt -71 

r e r e „e 

! - dz = | - ; - dx 

J . z J x + 17 

“00+ I7 —00 


izt 


and for convergence we must assume that t > 0 . Considering J - dz round 

a contour like that in Figure 8.3.1, except that the segment [-R + i 7 , 

R + i 7 ] is in the upper half-plane, 


, 00+i T izt 

1 r 0 


ETT J 


dz = 0 . 


•00+17 


If t < 0, then we consider the semi-circle in the lower half-plane. In this 
case the integrand has a simple pole at z = 0 within the contour and hence 
one shows that in this case 


. 00+i7 izt 

1 r e 


iTT I 


2n1 J . Z 

—00+ i 7 


dz = -1. 


When 7 < 0 we prove similarly that 

00 +i T izt 1 (t > 0 ) 

A j - dz = { 

Zn 1 J 


—00+17 


0 (t < 0 ) 


Interpretation 

If G(z) = -r~ j2n * z t ^ ien has, from the point of view considered in §8.3, two 

domains of analyticity, viz. the upper half-plane and the lower half-plane. 
When G(z) is taken together with the upper half-plane one would expect it to 
be the Fourier transform of 
0 (t > 0 ) 

f(t) = 

-1 (t < 0 ) 

and it is easy to check that this is indeed the case. 


When G(z) is taken together with the lower half-plane one would expect it to 
be the Fourier transform of 
1 (t > 0 ) 

f(t) = 

0 (t < 0 ) 

and again it is easy to check that this is the case. 

8 . If a < 7 < 7 2 < b we consider for R > 0, J G(z)e 1Zt dz round the 
rectangle with vertices - R + i 7l , R + i Ti , R + -R + i 7 Then 
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R+i 7 ^ 

|J G(z)e^ Z ^dz 
R+i 7l 


N <l'' 1 l + lT 2 l> - |t| 

< g.e Ai 2 - 7 X > 


-> 0 (R oo) ; 

and a similar result holds for the integral over [-R + i 7 ^, -R + i-y,,] 
Since the integrand is analytic within the rectangle we see that 


R+ 1 7 


R+ i 7 


lim 

R^-oo 




. izt, 
G(z)e dz 


— R+ i 7 


lim 

R->oo 


J 


„, v izt, 

G(z)e dz. 


-R i 7 


9. One checks that G(z) is properly defined for Im z < b by noting that the 
integral converges. To prove analyticity we use the result of Theorem 
4.2.3. 


[Remark: The convergence of Definition 4.2.2 is usually described as local 
uniform convergence.] 

Let K be a compact subset of {z: Im z < b], Then there is some v > 0 so 
that b - Im z > r; for all z e K. Set 


izt 

f (z) = f f(t)e dt 
n J 

0 

and then for m > n and z e K, 

m -izt 

|f (z) - f(z) I = IJ f(t)e dt| 
m n J 

n 


< K J e bt .e yt dt (y = Imz) 

n 

m _ t- 

< K J e V 'dt 

n 

K -r? n 

< - • e 

v 

K ~ nn 0 

Given e > 0 choose n„ so that - e < e and then 

0 v 

|f (z) - f (z)| < £ (z €K f m f n > n ), 
m n 0 

and so the sequence (f ) converges uniformly on K. Hence the result follows 
from Theorem 4.2.3. 

10. Note that here and later u(t) is the unit step function, also known as 
the Heaviside function, defined in Example8.4.1 at the bottom of p.373. 

Note also that the next few problems are meant to illustrate the way in 


'S 


$ 
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which the theorems of §8.4 may be used. However, even if one or other of 
these theorems is applicable it may be easier to argue directly. 

We have that 

00 

F[u(t).t ] = J t e dt. 

If z = x + iy, 

■ -izt. | -ixt+yt, yt 
I e I = I e I = e-* , 

and so for convergence of the integral we require Im z = y < 0, 


1 0 

Suppose that z = re with r > 0 and -n <6 <0, so that Im z < 0. Consider 

cx — i zc 

i e integrated round the contour consisting of the real axis from 0 to 

R > 0, the arc of the circle {|? j = R} from ? = R to { = Re"^ 0 7 r /2j (note 

-n <6 < 0) and the segment of the ray arg ? = -0 - n - from Re^ 6 7r / 2 j to 

the origin. If a a > 0 is not an integer we would indent the contour at the 

origin. ? 01 = e a ^°^ i s determined to be real on the positive real f-axis 
and elsewhere within and on the contour by continuity. 


ex — i zf 

Observing that ? e has no singularities within the contour and using the 

calculus of residues and letting R -* oo we find that 

00 00 

r a -izt-^ r a ia{-8-n/2) -rt l (-0 -tt /2) _, 

| t e. dt = J t e .e .e dt 


0 


0 


= e i(o + l)(-»-T/2) J t“e' rt dt 

0 


i(a+l)(— 8— n/2) 1 (• a — U 

= e .-- | u e du. 

a + 1 J 
r 0 


, . i<9 a + 1 

(ire ) 


r (a+1) 


r(a+l) 


( iz) 


a+1 


We have used that as R ^ oo the integral along the arc of {|? | = R} tends to 
0 and that we can let the indentation, if one is necessary, shrink to 0. 
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11. We consider only 


- i. z t 

F[u(t) COSh u t ] = — r~— f COSh wt.€ dt 

^ 2n 0 


00 ut 


- 727 5 


e + e 


-izt,, 
-.e dt. 


Now, assuming u is real, 

00 

(w-lZ)t,. 1 


J e 
0 


dt = 


u - 1Z 


J e -'" +iz)t dt = 

0 


u + 1Z 

with Im z < - |w|. Hence 


F[u(t) cosh ut] = 


r-> 


2-J2tt lj - iz u> + iz 


1 iz 

J2n 2 


2 * 


LJ + Z 

Remark: If we observe that cosh ut = cos iwt we can use the result of 
Example 8.4.5 on p.374. However, one has to be careful about the associated 
strip of analyticity. 

12. Here we can proceed as in the preceding problem. 

13. The transform asked for is 


w # yo 

757 J te dt =7^J ‘ 


(a-iz)t 
e dt. 


With the appropriate restriction on Im z the integral is easily evaluated 
using integration by parts. 

14. One can deal with this problem in the same way as the preceding one. 

1 


15. The function 


z ( z 2 - 1) ( z 2 + 1) 


has 5 singularities and these are at 


z = 0, ±1, ±i. Considering the discussion in §8.3 if we consider this 
function as a Fourier transform there is an associated strip of analyticity 
which should emerge from the solution. 


In partial fractions 

_ 1 _ 

z(z 2 - 1) (z 2 + 1) 


1 1 _ 

z + 2 2 


z - 1 
1 


2 2 
z +1 


1 1 , 

- + - ( 
z 4 1 + z 1-z 


1 ) + 


4i 1 - iz 1 + iz 


If ? = (, + itj , then 
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? t-izt 


-l r s u —1 r (f — i z ) t 

757 J e dt = 757 I e dt 


0 

1 1 
\J2n *? - iz f 

provided £ + y < 0, i.e. Im z < - £ and so 

ft 1 1 

F[ u (t) e j = - -f—. - ±-r- 

\12 n ( - 1Z 

provided Im z < - Re { . Hence 

f(t) = -i72» u(t) + j i72ir u (t) e” 1 * - | (- i72, u (t) e lfc ) 

+ -Jr- (- 72ir u(t) e*) - -ji (72, u(t) e t ) 

e _it e U It 1 -t 

= i' 72 ' u <t> [- 1 + —r + — - 4 e + j e k ] 

= -i\/27r u(t) [1 - | cos t + | sinh t]. 


16. From the discussion at the beginning of §8.5 one finds that 
F[y] = 


F[e '* ] 


(iz) + 3(iz) + 2 

_ F[e~^ ] 

(1 + iz )(2 + iz)’ 

From solution to P.5, 

_ r — | t | _ 1 , 1 1 x 

F[e 1 = 727 ( “ + -- rz) 


IZ 


IZ 


and so 


F[yl " 757 ( I“7‘iz + l Ti 1 ( I /iz ” 7+ iz 


nr) 


■J 2n 


[• 


+ \ ( 


1 + 


■J 2n 


[ 


(1 + iz) 

- ( 

_JL_ 

(1 + iz) 2 


2 1 + iz 1 - iz 


1 + iz 2 + iz 
1 1 


nr) " o ( 


1 


1 1 
+ -. 


3 1 - iz 2 + iz 

1 


nr)] 


2 

+ 


2’1 + iz 6*1 - iz 3*2 + iz 




Hence 


y(t) = [(t - |)e 1 + | e 2t ] u(t) + | e 1 [1 - u(t)]. 

It is easily checked that this is a solution of the equation for t > 0 and 
for t < 0 . 

17. From Example 8.3.1, p.366, 

1 w 


F[u(t) sin wt] = 


J2 n 2 2 

U) — Z 


(Im z < 0 ). 
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Hence 


r n 1 W 1 

i!yJ = - 727-T- 2'—2 ---"• 

w - z z - 3iz - 2 

One can now proceed as in the preceding solution or one can argue directly 
and we shall do the latter. 


The relevant strip is the lower half-plane so we consider 

1 w 


, oo+i 7 izt 

1 r 1 w e 


Y () 727 J* . 72 rr ‘ 2 2*2 

—00+17 w — z z — 3 iz — 2 


dz 


00+ i-y 


2n 


J 


izt 


2 2 2 

•00+i 7 ( w — z )(z - 3 iz- 2 ) 


dz, 


where 7 < 0. Using the contour of Figure 8.3.1 on p.366 and arguing in the 
usual way we find that for t > 0 , 


00+ iy 


izt 


J 


2 2 2 

•00+17 (w - z )(z - 3 iz - 2 ) 


dz = 2?r i x sum of residues at ±w , i, and 2i 


iwt 


-iw t 


-t 


■21 


= 2rri [ 


(— 2w ) (w --3iw—2) 2w (w +3 iw-2 ) (w +1)(—i) (w +4)(i) 


n 1 e 

77 ~2 


•iw t 


iw t 


+ 3 iw —2 


-3iw—2 


■] + 2n [ 


-t -21 
e e_ 


2 , 
w +1 


w 


' + 4 


2 ' 2 
+ 1 w +4 


. , 2 -iwt iwt v , -iwt iwt, -t 

rr 1 (w -2 ) ( e -e ) -3 iw ( e +e ) . , „ , e , e 

- [- X - 55 - J + 2 JT [--- + 

(w -2) +9w 1 

2 -t -21 

n 1 -2i(w -2)sinwt-6lwcoswt. _ r e e 

— [ - r~r - 1 + 2 ’ ft + ~i 


-21 


+ 5w +4 


w +1 w +4 


2 -t 

2ir . (w -2 ) sinw t + 3 wcosw t, _ r e 

“ [ -TTT-- 1 + f- — 


-21 


■]. 


w +5w +4 
Hence for t > 0, 

^-t 

y (t) = w [ 


2 2 
w +1 w +4 


w +1 w +4 


-21 2 

e n r (w - 2 )sinwt+ 3 wcoswt n 

j [ ~ ~ J • 


4 _ 2 
1 +5w +4 


The general solution of the homogeneous equation is 

-t -21 

y(t) = Ae + Be 

and so the general solution of the given equation for t > 0 is 


-t - 2 t r (w - 2 )sinwt+ 3 wcoswt, 

y(t) = Ae + Be - [± -- - -]. 

w +5w +4 
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[Remark: Had we noted the form of the general solution of the homogeneous 

■> ■ “t -21 . 

equation earlier we could have dropped the terms involving e , e .J 


To find the solution of the problem we determine A, B from 

y (t) 0 (t -> 0 +), y'(t) -► 0 (t 0 +) and hence 

2 2 
w (w +4) _ w (w +1) 

A = —:-— / B — — 


4 _ 2 . 
> +5w +4 


4 _ 2 . 
w +5w +4 


18. If we apply the Fourier transform to both sides of the equation we get 

[ ( i z) + . . ~~r ] F[y] = F[u(t) e t ], 

(iz) 

where the second term on the left is the result of an appropriate 
integration by parts. Hence 
1 iz 


F[y} = 


727 


( 1 +iz)( 1 -z ) 

1 

727 1 2 * iz+1 4 * iz+i T 4 * iz-i 


1 r. i 1 + + 4i.r=4] 


and so 


1 -t 1 +i -it 1 -i it 

y(t) = "2 e + — e + -4“ e 


= lt _e + cos + sin ^1* 


Considering the homogeneous equation we expect that it can be differentiated 
to give 


at 2 


+ y(t) = o 


which has solutions of the form A cos t + B sin t. Considering 
y(t) = cos t, y(t) = sin t in 

^ + J y(r )dr = 0 
dt 0 

we see that a solution is of the form A cos t and so we are looking for a 
solution of the original problem of the form 

y(t) = - | e | sin t + A cos t. 

If y (t) -* y (t + 0) , then 

y 0 = - 5 + A, ft = y Q + j 
and a solution of the problem is 
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y(t) = - 


1 -t 1 . ^ l x 

- e + - sin t + (y n + -) cos t 


The next 3 problems seem rather difficult and are not comparable in many 
ways to the illustrative examples of §8.6. We shall, therefore, be 
concerned in the following solutions with a general discussion of the 
problems rather than with finding explicit solutions in all cases. 

19 . Following the argument of Example 8 . 6.5 we are solving the equation 


2 2 
<9 u 1 d u 

2 ~ 2 2 

d X a d t 


- bx(n- - x) (0 < x < 7T , t> 0) 


subject to the conditions 

u(x, 0) = 0, u fc (x, 0) = 0; 

u(0, t) = 0, u(?r , t) = 0. 
If we consider Im z < 0 and 


1 ~ i z t 

U = U( x, z) = F [ u ] = J u(x, t) e dt, 


z b 

U + — U = - -rr-r- x Or - X . 

xx a \/ 2 n . iz 

The general solution of the homogeneous equation is 

Z X Z X 

U(x,z) = A sin — + B cos ——, 
a a 

where A, B are functions of z. A particular solution of the above equation 
2 

is Cg + c^x + c 2 x with 


2C 2 + h <C 0 + C 1 X + °2 x2) = ' 

a 


X (rr - X ) 


so that 


2 

z 

b 


a 2 b 

~2 C 2 
a 

J2n . iz ’ 

C 2 

3 ' 

•J2n . i z 

2 

z 

b TT 

pi 

a 2 b^ 

~2 °1 
a 

J2n . iz ' 

c i 

3 ' 

J2n . iz 

2C 2 + 

h c o = °' 

°o 

2 a 2 

■ - 2 c 2 
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Hence the solution is 


Z X zx 

U(x, z) = A( z).sin — + B(z) cos — 

a 2 b 


\/2n . iz 3 


,2a , 2 

(—- + ;rx - x ) . 


From the boundary conditions we can determine A and B and then use the 
inversion formula. 

20. This problem involves solving 


2 

d U 


d X 


, 2 
1 d u 

~2 2 
a <91 


(t > 0 ) 


subject to the conditions 

u(x, 0 ) = 0 , u fc (x, °) = °' u (°' = 9(t). 

Because of the first two conditions one might consider taking a one-sided 
Fourier transform with respect to t. However, if one does this one arrives 
at an impasse because one is unable to determine a function A(z) in the 
general solution by the given conditions. In fact it seems necessary to use 
the Fourier sine transform as in Example 8.6.4 on p.384. Consider 


U = U(z, t) = - V- i f u(x, t) sin zx dx. 
* 0 


Then, as on p.385, 
U 


tt 


2 2 ,2 . 2 ... 

+ a z U = - v/-. iza g(t) 


and so, using the result on p.385, 

t 2 2 

,2 . 2 , , . -a z (t-r ) 

U(z, t) = A sin azt + B cos azt - V- iza J g(r ) e Qt • 

" 0 

The condition u(x, 0) = 0 requires that B = 0, and the conditions 

u(x, 0) = 0, u(0, t) = 0 imply that g(0) = g(t) = 0 so that = Aaz 

when t = 0 and therefore A = 0. Therefore as on pp.385, 386 we find that 

2 


-] 


g(r ) exp[—-— 
x r 4a (t-r) ^ 

u(x, t) = -z—j- J -3/2 dr • 

2aJn 0 (t-r) V2 

[Remark: Verification of this solution is really rather tricky.] 

21. If one follows the hint then one is concerned with a doubly infinite 
string whose displacement satisfies the wave equation 


2 

d u 
2 

d X 


, 2 

1 d_u_ 

~2 .2 * 

a d t 
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This means that the extended problem must be one in which the string is 
freely vibrating for t > 0. This assumption is physically plausible given 
the hypotheses of the problem. 

We consider Example 8.6.1 with f(x) of the example replaced by f(x) (x > 0), 
-f(-x) (x < 0) and g(x) = 0. The solution then, from p.382, is 

u(x, t) = ^[f(x + at ) + f ( x “ at )]' 

where is the function of the problem extended to an odd function. 

For the final part of the problem consider, as in Example 8.6.4, 
u(x, t) = u 1 (x, t) + u 2 (x, t). 

Since f(0) denotes a displacement of the whole string up or down we can take 
£( 0 ) = 0 and consider 

u (x, 0 )=f(x), ~ u (x, 0 ) = 0 , u ( 0 , t) = 0 (i.e. the problem above); 

1 a t 1 1 

u_(x, 0) = 0, ~ u (x, 0) = 0, u (0, t) = g(t) (i.e. P.20). 

2 a t 2 2 
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Chapter 9 

Laplace Transforms 


| READ §9.1 | 


The Laplace transform is regarded as a two-sided transform of functions 
which are 0 for t < 0. Whilst this means that one can use the results of 
Fourier transforms rather directly it might have been better to give an 
account that restricted attention to the one-sided case as this is an 
essential feature of the Laplace transform. 

Definition 9.1.1 The transform H(z) = L(f) will be analytic for Re z > b. 

Theorem 9.1.1 A Laplace transform H(z) is always analytic for all z of 
large enough real part and so the inversion formula will be valid for all 
large enough positive 7 . 

Self-Assessment Questions 

SAQ 1, 2, 3 are Exercises 9.1, #1, 2, 3. 


| READ §9.2 | 


One should not be too concerned about the proofs of the theorems in this 
section. If one proceeds in a formal manner one arrives at the conclusions 
of these theorems. In Theorem 9.2.2, for example, we have, assuming that 
the symbols have a meaningful interpretation, 


L(f) = f £(t) e Zt dt = - f(t>e 

0 z 


-zt 


+ - f f'(t)e Z ^dt 
z ** 

0 0 


= + i L[f<t)], 

z z 

and this is the formula of the theorem. In any special case it should be 
clear whether the steps above are valid and the result true in the special 
case. 


Example 9.2.8 The notation used in this example is not in conformity with 
more or less universal practice and so is very misleading. If one has a 

derivative f^ n ^(t) then f^ n ^(kt), as one would expect, is the value of the 
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nth derivative at kt, i.e. we differentiate n times then give the variable 
the value kt. We do NOT give the variable of f the value kt and then 
differentiate n times with respect to t. In other words 


f^(kt) = —f (x) | and NOT —— f (kt). 

n - n 

dx x=kt dt 

However, Dettman uses the latter notation so that 
t J Q "(kt) + J 0 '(kt) + k 3 t J Q (kt) = 0 
does NOT mean what one expects, but 

72 J o (kt) + a¥ V kt) + k2fc J o (kt) = °- 

dt 

Perhaps it is also worth noting that one could take k = 1 initially and then 
find L((kt)] by a transformation of the variable in L[J (t)]. Note also 

that Jg(-t) = J Q (t) and therefore we can assume that k > 0 without loss of 

generality. 

Self-Assessment Question 
SAQ 4. is Exercises 9.2, #1. 


| READ §9.3 | 


Theorem 9.3.1 In the displayed formula of the proof the minus signs should 
be plus signs. 

Theorem 9.3.2 et seq. One should read these theorems in order to get an 
idea of how one deals with the inversion problem. It is normally better, 
provided one has a grasp of the general ideas of these theorems, to treat a 
particular problem on its own merits. 

Problems for §9.1, §9.2, §9.3 

PI, 2, 3 are Exercises §9.2, #2, 3, 4. 

P4, 5, 6 are Exercises §9.3, #1, 3, 4. 


| READ §9.4 | 
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Note carefully the general method for using Laplace transforms in order to 
solve ordinary differential equations. The examples are rather complicated 
and you should not spend too much time on the details unless you are 
interested in the physics context. Observe how the Laplace transform copes 
with initial conditions. 

Problems for §9.4 

P7, 8 , 9, 10 are Exercises 9.4, #1, 2, 3, 4. 


§9.5 is a special topic and does not form part of M828 . 

% 


| READ §9.6 | 


Note that this section does not contain anything that is essentially new. 
As the introduction says, this section gives the Laplace transform version 
of what was done in § 8.6 using the Fourier transform. 

Problems for §9.6 

Pll, 12, 13, 14 are Exercises 9.6, #1, 2, 3, 4. 


| READ §9.7 | 


Problems for §9.7 

P15, 16, 17 are Exercises 9.7, #1, 2, 3. 


Solutions to Self-Assessment Questions 
1. We have 

00 00 

L[u(t) ] = f e Zt dt = - - | = - (Re z > 0). 

J „ z „ z 

0 0 

The verification of the inversion formula requires us to think of L[u(t)] as 
a two-sided transform and so we are concerned with the formula at the bottom 
of p.399. Hence the inverse transform is 

7 + ioo f t 

^—r J d? (Re 7 > 0 ) . 

Zn 1 •> , f 

7 = loo 

Iff = £ + iv, then 
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(t 


. i it 

e j = e 


1 ')' + i R 

and so for t > 0 we consider —-r f together with the integral over the 

2 n 1 J _ 

7~lR 

semi-circle from 7 + iR to 7 - iR on the left. As R -> 00 , the integral over 
the semi-circle tends to 0 and so 

. ft 

_ f ~= Rps - 

2 ? 


, 7 + ioo f t 

1 r e _ _ e 

—— 1 - df = Res — 

In 1 J . f 5 

7 — loo 


? =0 


If t < 0 we use a semi-circle on the right and note that in the contour the 
integrand has no singularities. Hence 

A- r 

2 tt i j 


7 + ioo f t 


* ^ j 

- df = 

f 

7 - loo 

0. 

, 1 

1 we have -—- 

2 rr 1 • 

7 + ioo 

r 

7 = ioo 


d? 


which does not converge. However, if one 


considers Theorem 8.2.1 and its proof one sees that it is the Cauchy 
principal, value that is relevant. Now 
7 + iR 


-K r 

2 n i J 


df 


f 2 n i 


- [log( 7 +iR) - log( 7 -iR)], 


7 -iR 

where we determine the logarithms to be real,for R = 0 and continuous in R. 
Then 

Im[ log (7 + iR) ] -* ?r/2(R -> 00), Im[ log (7 -iR) ] * ~n/ 2 {R -► <*>) 

and 

17 + iR 


log |~— 7 -pl * 0 (R 00 ) 

7 — IK 


and so 


7 + iR . , 

f 3L * 1 [“ - <- il’] = - 

!!;r i 1 . f 2 n i 2 2 2 

7 — iR 

u( 0 + ) + u (0 ) 


2. Assuming w is real we have 

00 


— Z t 

L[u(t) cos ut] = J cos wt.e dt (Re z > 0 ) 

0 

00 iu t — iw t . 

r e + e -zt 

- j -5- e dt 

0 


„ (iw-z)t -(iw+z)t 

- - - ® — - ] | 

2 lw-Z 1W+Z 


0 
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= - [-r-i- . ■■■ — ] 

2 iw + z i(jj—z 


2 2 * 
z +w 


For t > 0 we consider the usual semi-circle on the left from 7 + iR to 7-iR 
and obtain, for 7 > 0, 

, 7 + ioo , Zt Zt 

1 r z zt, „ ze 1 , n ze 1 

— J — - e dz es — - I . Re 2 2 ' 

7—loo Z +w Z +cj Z+lw Z +w Z = —loo 

, lw t , — 1 U) t 

iwe iwe 


lGJ t — lCJ t 

e + e 


= cos ut. 

The cases t < 0, t = 0 are dealt with in the standard way. 

3. Assuming that differentiation under the integral sign is permissible 


if LU(t>] - - at J fit) e ~ zta t 


= - J f(t) ~ [e Zt ].dt 

0 az 


= J t f(t) e Zt dt 
0 


= L[t f(t)3. 


From Example 9.1.1, 


L(u(t) sin ut] = 


2 2 
Z +w 


and hence 


L[t u(t) sin wt] 


— [—-—] 
dz 1 2 2 J 

z +w 


2 2 2 * 

(z +w ) 

4.(a) Assuming a is real, 

0 ° 

L[u(t) t n e at ] = J t n e at e Zt dt (Re z > a) 

0 


= j tV (z " a)t dt. 

0 

Now suppose that z = x is real and change the variable of integration by 
(x - a)t = t and then 
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L[u(t) t n e at ] 


1 


r n -r 
r e dr 


(x-a) n+1 0 


n! 


(x-a) 


n+1 


n! 


By analytic continuation the transform is ^ 

(z-a) n 

(b) Assuming w is real, 


(Re z > a). 


— 21 

L[u(t) cos h wt] = J cos h wt e dt 

0 


0 ° 

1 p . w t o t 2t . I j v 

= - j (e + e ).e dt (Re z > |w|) 
A J 


o 


, (w~z)t ~(w+Z)t 00 

= i [S_2_] | 

2 u-z w + z q 


- [ 1 
^ L 


2 w+z u—z 
z 


2 2 
Z —w 


(c) Assuming b is real and a > 0, 
bt °° 

L[t u(t~a)e ] = J t u(t-a)e 

0 


-(z-b) t 


= J te 

a 


- (z-b) t 


dt (Re z > b) 

-(z-b)t 


I + -— f e 

z-b J 

a a 


dt 


t -(z-b)t 


z-b 


■ e | + -itr f e "'’dt 

z-b J 


-(z-b)t 


z-b 
a a 


a -(z-b)a 
e 


z-b 


1 -(z-b)t 

-r- e 


(z-b) 


a -(z-b)a , 1 -(z-b)a 

• e + - — e 


z-b 


(z-b) 


ab -az ab -az 

e .e ae .e 

--- -f 


(z-b) 


z-b 


—7 — t 

(d) The function J e cos w(t - r )dr is the convolution of e and cos wt, 

0 

The one theorem it is perhaps useful to keep in mind from §9.2 is Theorem 
9.2.7 which deals with convolutions. From this theorem the transform asked 
for is 
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00 


00 


— f — ( 2 4- T \ fc “Zt 

L[u(t) e ].L[u(t) cos ut] = J e dt.J cos w t.e dt 

0 0 


1 z 


z+1* 2 2* 

Z +w 


[The answer given by Dettman is wrong.] 


Solutions to Problems 
1. For x > b, S >0 consider 


(x) = J dt J f(u) 


-tu 


du. 


The hypotheses of the problem ensure that I (x) exists for x > b 

6 


and it can be shown that for such x, 6 the orders of integration 
interchanged. Hence 


and S > 0 
can be 


-tu 


I (x) = J f(u)du J e dt 


= f £i-Sl.e' XU du. 

s u 


The condition on 


f(t) 


as t -* 0 shows that 


lim 


V s0 “ «“o V x> - J u 


£<U) .e- XU du = • 


In the integral for I (x), given the hypotheses of the problem, one can 

0 

prove that can be interchanged with the integration with respect to t 


and hence 


I Q (x) = J G (t)dt = H(x). 
x 

Consequently H(z) and L[ ^^ ] as a function of z are equal for each real 

z > b and so the result follows by analytic continuation. 

2.(a) Let w be real and put 

OO 

H(z) = f Sin e zt dt (Re z > 0). 

0 k 

Assuming that differentiation under the integral sign is justified (which it 
is) , 
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H'(z) = - J sin wt.e Zt dt (Re z > 0) 
0 


2 2 
Z +w 


Considering z = x > 0 and noting that H(oo) = ■ L: *' ITl H(x) = 0, 

x-»-oo 


H( x ) = J H' (t ) dt = J 


,2 2 

x t +w 


dt 


“1 LJ 

= - tan (-). 

x 


By continuation 


H(z) = tan ^(-) (Re z > 0). 
z 

(b) The proof is similar to that of (a). 

(c) Let 

°° J l (t) -zt 

H(z) = f —-—.e dt (Re z > 0). 

o 1 

For convergence at the upper limit consider the representation for J^(t) on 

p.316 and for convergence at the lower limit consider the series for J (t) 

on p.312. These results can also be shown to justify differentiation under 
the integral sign and so 


H 


-zt 


' (z) = - J J (t)e dt = - [1- - 

0 y/l + Z 


]» 


from the result at the top of p.410. 

Taking x > 0 and integrating from x to oo and noting that H(R) -> 0 (R <x>, R 
real) we obtain 

R 

I 

R-^oo 


- H( x ) = ~ J [1- 

X y/l + t 2 


1 dt 


lim r „ , . 2 , 2 v , 

[R - x - (Vl + R - \/l + x )]. 

R-*oo 


Since 


lim r „ . 2. lim r 

R*oo 


R + y/l + R“ 


■] = 0 , 


we deduce that 


H(x) = \/l + x' 


x. 
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In the integral for H(x) put t = ku with k > 0 and then 


°° J i (ku) -vku 
H(x) = J -A-.e XKU .k du 


0 


ku 


oo J (ku) 

r- 1 “XkU.. 

-.e du. 


0 


u 


From this result and analytic continuation, for Re z > 0, 

J l (Kt) 99 

L[— E -] = H(|) = ± [^(k^ + z) - zl. 

Since J^-kt) = - J^fkt) this result is also valid for k < 0 and letting 

k -> 0 one checks the result for k = 0. 


Note that \/(k + z ) in {Re z > 0} is positive on the positive real axis and 

is determined elsewhere by continuity. 

3. The reference to Theorem 9.1.1 seems to be a mistake. We suppose that 
L[f] exists and that Re z > 0. Let N be a positive integer and then 

L[f] = n'i r f(t) e_zt at 


.. N—1 (k+1 )t 

llITl V f ~ Zt AX- 

= „ ) f(t)e dt. 

N->oo L J 

k = 0 KT 

In the kth term of the sum put t = kr + u and note that by periodicity 
f (kr + u) = f(u). Hence 

(k+1 )r 


~ ZKr r C, X ~ ZU A 
f(t)e dt = e J f(u)e du 

0 


— zk; 


kr 


and so 


Nr N-l r 

J f(t)e Z dt = [ ^ e Zr ]J f(u)e ZU du 
0 k = 0 0 

, -ZrN r 

1 - e r s ~ zu a 

= -. f(u)e du. 

, “Zr J _ 

1 - e 0 


-Zr N 

Since Re z > 0 , e ■* 0 (N -» oo) and therefore, letting N oo, 


. t 

L[f] = --- f f(u)e ZU du. 

, —Zr J _ 

1 - e 0 

One would expect the rider to be a special case of the general result, but 
it is not. What would seem a more natural choice is the fractional part of 
t, i.e. {t} = t - [t]. However, for the question asked, 

[t] =n(n< t < n + 1) 
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and so 


00 


J 

0 


[t]e zt at 


oo n+1 

J n | e Z dt (Re z > 0) 
n=0 n 


r n r -zn -z(n+1), 
l -le - e ] 

n=0 


, — Z oo 

1 - e r -zn 
- \ ne 


n=0 


-z 


1 - e d . r -zn, 

—i—• ai ‘ L e 1 

n=0 



1 

z(l - e _z ) 

4. This problem seems unduly elaborate for an exercise that is meant to 
illustrate the general ideas and results. One must find the partial 
fraction form of the given function. Since 


z 2 + 2z + 5 = (z + l) 2 + 2 2 = (z + 1 + 2i) (z + 1 - 2i) 
we consider 


B 


D 


E 


(z+1)(z—2) 2 (z 2 +2z+5) 


z+1 z-2 

Multiply by z+1 and put z=-l and then 


(z-2) 


z+l+2i z+l-2i* 


(— 1 — 2 ) 2 [(— 1) 2 + 2 .(— 1 ) + 5 ] 

-1 

9(1^2+5) 

1 

36* 

2 

One can determine D, E similarly. Multiplying by (z - 2) and putting z = 2 
determines C. Finally B can be determined by multiplying by z and letting 
z ■> oo to obtain 

0=A+B+D+E, B=-A—D— E. 

The relevant function can then be found as in Example 9.3.1 on p.412. 

5. We shall verify entry 13 only. Because §9.3 deals with inversion in the 
manner it does and G(z) is on the left of the table you might think that you 
are required to start with G(z) and find f(t). However, in a verification 
one should proceed in the simplest possible way and so in this case one 
should start with f(t) and find G(z). 
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If 


then 


r 1 - COS wt -Zt,, 

H(z) = f ---.e dt (Re z> 0), 

o 1 


H'(z) = - J (1 - cos wt)e Zt dt 
0 


.1 z . 

~ ~ z + ~2 2 * 

z +w 


Considering z = x > 0 and noting that H(R) 0 (R ->■ <x>) we get 


lim 


- H(X) = J H ' (U) dU 


= - i [log(x 2 + u 2 ) - log x 2 ] 

2 2 

lim 1 . R +uj 
+ R*» 2 109 


and so 

H(x) = | [log(x 2 + w 2 ) - log x 2 ]. 

The result now follows by analytic continuation., 

6 . Note that b 2 should be • A cubic polynomial can be put in 

the form 

ftz{(z - a) 2 + b 2 } + /?{(z - a) 2 + b 2 } + 7 z + 5 

p ( z ) 

and hence one obtains the given form of ——- with A„ = 0. 


Q(z) 


To find f(t) so that 
P(z) 


L[f(t)] = 


Q(z) 


we use 8 and 9 of the table on p.415 and the results obtained from 8 and 9 
by differentiating with respect to z. 

7. If Y(z) = L[y(t)], then 

r>o 

(z+1)t 


J y"(t)e Zt dt + 3 J y'(t)e Zt dt + 2 Y(z) = J e' 
0 0 0 


dt. 


Now 


i X . 

J y'(t)e Zt dt = y(t)e Z | + z J y(t)e Z dt 

0 0 0 

= - y Q + z Y(z); 
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00 oo 


J y"(t)e Z ^dt = y'(t)e \ + z J y'(t)e Zt dt 

0 0 0 


y 0 " Y o z + z Y(z) * 


Hence 


z 2 Y(z) - y Q z - y^ + 3z Y(z) - 3 y Q + 2 Y(z) = 


and so 


Y(z) = 


1 /z+l + y Q z + Yq + 3y Q 


z + 3z + 2 
1 


v z + v' + 3 v 
y 0 y 0 y 0 


(z+1)(z+1)(z+2) (z+l)(z+2) * 

We now determine y(t) from the inverse Laplace transform. 
8 . If Y(z) = L[y(t)] then 

00 00 

J y"(t)e Zt dt + w Y(z) = J sin wt.e 2 dt. 

0 ° 0 

Now 


00 OO 


J y"(t)e zt dt = y'(t)e zt | + z f y'(t) e zt dt 

0 0 0 


-z t i 2 

= z y(t)e | +z Y(z) 


= z Y ( z ) ; 


and 


r • . ~zt w 

sm wt.e dt = —-- 

0 z +w 


Hence 


Y(z) = 


. 2 2 2 2 
( Z +w ) ( Z +w ) 


For w ^ w^ we use 4 of the table on p.415 and when w = w^ we consider the 

differentiated forms of 4 and 5. 

9. Let Y(z) = L[y(t)] and consider 

oo oot oo _ . , 

J y'(t)e Zt dt + J [J y(r )dr ] e 2 dt = J e Z+ dt. 

0 0 0 0 
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Now 


f y'(t)e Zt dt = y(t)e Zt | + z Y(z) 

0 0 

= - y o + z Y(z); 

°0 t _ , t oo 

J [J y(r )dr ] e Z dt = - - [J y(r)dr].e | + - Y(z) 

0 0 z 0 0 

■ ; Y(z) - 

Hence 

- y Q + z Y(z) + \ Y(z) = ^ 
and so 

z(z+l+y ) 

Y(z) = ---. 

z +1 

One now use the inverse Laplace transform to find y(t). 

10. Let X(z) = L[x(t)]f Y(z) = L[y(t)] and consider 


Now 


00 _ 00 _ , 

J x'(t)e dt + J y '(t )e Z dt - X(z) + 3Y(z) = 

0 0 

00 00 _ 

f x'(t)e" Zt dt + J y'(t)e Z dt + 2 X(z) + Y(z) = —j 
0 0 


f x'(t)e Zt dt = x(t)e | + z X(z) 

0 0 

= - x 0 + z X(z); 


-zt 


and 


-zt 


f y(t)e dt = - y + z Y(z) 
0 


Hence 


(z - 1) X(z) + (z + 3) Y(z) = 


FTT + x o + y o 


(z + 2) X(z) + (z + 1) Y(z) = T - rJ + x Q + y 0 . 

One now solves for X(z) and Y(z) and then one uses the inverse transform to 

obtain x(t) and y(t). 

11. The complementary function is 

zx/a _ -zx/a 
Ae + Be , 
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where, remember, A and B are functions of z. Denoting — by D we have 


(D 2 - ~)U --1 g(x) 

a a 

and so a particular integral is 
z 1 


2 2 2.2 
a D -z /a 

= - -i t 1 


• g(x) 


2a D-z/a D+z/a 


] g(x) 


2 a 


1 . zx/a r -zt/a 

[e j e g (t)dt - e 

0 


-zx/a j. e zt/a g(t)dt]t 
0 


The general solution, therefore, is 


X X 

zx/a „ -zx/a 1, zx/a r -zt/a -zx/a r zt/a 

U(x,z)=Ae / +Be ' - —[e j e g(t)dt - e J e g(t)dtj. 


From U(0, z) = 0 we see that 
A + B = 0. 

Since 

x 


1 


u (x, 2 ) = [A - — J 


-zt/a 


g(t)dt] e 


zx/a 


+ [....] e 


0 


-zx/a 


if U(x, z) 0 (x oo) for Re z > c > 0, then 


00 

A(z) = / J e zt a g ( t ) dt . 

2a 0 

[Note that there is a misprint and x -* 0 should be x -* oo. ] . Hence A, B are 
determined and the problem solved. 

[Remark: This solution agrees with that given in Example 9.6.2 which is the 
same problem.] 

12. The conditions are those of Example 9.6.2 with u(x, 0) = g(x) replaced 
by u(x, 0) = 0 and u (x, 0) = 0 replaced by u (x, 0) = h(x). In this case 

if 


OO 

U (x , z ) = J u (x, 
0 


t)e Zl: dt 


then 


U 

xx 



h(x) 

2 ' 


a a 

Hence the solution will be the same as on p.430 except that g(x) will be 
replaced by and so, noting that we put h(t) = -h(t) (t > 0), 
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U (x, z) = i.- f h(x + ar/ )e Zr? dr/ + -. - f h(x - at?)e Zr? djj. 

2 Z 0 2 Z 0 

Now L[1] = - and hence from the convolution result we obtain 
z 

t t 

1 1 

u(x, t) = - J h(x + ar)dr + - J h(x - ar)dr. 

0 0 

13. Since the wave equation is linear if we can find solutions u^, u^, u^ 
satisfying appropriate initial and boundary conditions so that 
U = U 1 + U 2 + U 3 

satisfies the conditions of the problem then this u will solve the problem. 
Consider then the following problems: 


(i) 

u 1 (x, 

0 ) 


u 1 (x, t) 

1 = o, u 

t = 0 

(0, t) 

= f(t) 

(Example 

9.6.1); 

(ii) 

u 2 (x, 

0 ) 

= g(x), 

d 

at u 2 (x ' 

o 

1 ! 

O 

II 

4 -> 

-P 

u 2 <0, 

t) = 0 

(Example 

9.6.2); 

(iii) 

u 3 (x, 

0 ) 


u,(x, t) 

1 = h(x), 

t = 0 

u 3 (0, 

t) = 0 

(Exercise 

9.6.2). 

14. We 

make 

the 

plausible assumption that if 

we consider a 

semi-inf inite 


string with x > 0 and initial displacement A sin rrx, then the part between 
x = 0 and x = 1 will vibrate as though it were fixed at these two points. 

We therefore have the problem of Example 9.6.2 with g(x) = A sin ?rx and 
since A sin n x is an odd function 

u(x, t) = | A sin [rr (x + at)] + | sin |>(x - at)] 

= A sin n x. cos an-1. 

This answer confirms that the above 'plausible' assumption is compatible 
with the mathematical model. 

15. By the method discussed at the beginning of §9.7 if U(z) = L[u(t)] r 
U ( z ) = L (1) + L [ t 2 ] U ( z ) . 

Here we have used that the transform of the convolution of two functions is 
the product of the transforms (Theorem 9.2.7). Since 

OO 00 

L( 1 ) = J e" Zt dt = L(t 2 ) = J t 2 e _Zt dt = — (Re z > 0 ), 

0 Z 0 z 

U(z) = \ U(z) 

Z 3 

Z 

U(z) = _JL_. 

z - 2 
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U(z) comes from an 'unknown* function u(t) and assuming it has a Laplace 
transform this will be analytic for Re z > c and some constant c. From the 

1/3 

form of U(z) one sees that one can take c = 2 . Note that a suitable 

value for c depends on the equation and not just on the half-planes of 
analyticity a rising from the particular functions that enter into it. 


Hence 


7 + ioo 2 

1 r z zt,, 

u(t) ■ ST J . -3—7 e dt ' 

7 — loo Z — Z 


1 yn ( 

where 7 is a number satisfying 7 > 2 . Considering - t- J together 

7 -iR 

with the integral over semi-circle on the left and letting R -* 00 we find 
that 

z 2 zt 

u(t) = Y[residues of —^-*e ]• 

z - 2 

1/3 1/3 1/3 2 

The residues arise from simple poles at z = 2 ,2 u, 2 u , where 

2n i/3 „ 

u> = e . Hence 

1/3 „l/3 -1/3 2 

l r 2 7 t 2 / wt 2 u t 

u (t) = ^ [ e + e + e ] 

r 2 k t 3k 

‘ Jo (3k)! ‘ 

[Remark: It is a useful exercise to verify this result.] 

16. In this problem u(t) is the 'unknown'. The right hand side of the 

equation is the convolution of u(t) and — and so 


+ iR 


F ( z ) = L[f (t) 3 = L[u(t)].L[ — ] 

t" 

00 —z t 


= U ( z ) . f --dt. 

J a 


0 t 

Taking z = x > 0 and putting r = xt we obtain 

°o —Xt . 00 . v-. 

f ®- dt = X _(1 a) J r ( a) e r dr 


0 t 


0 


= X (1 a) .r(1-a). 


By continuation we note that 
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.r(l-a) (Re z > 0) 


L( —) = z (1 a) 


where z ^ is taken to be real and positive on the positive real axis and 

is determined elsewhere in the right half plane by continuity. Hence 
1-a. 


U(z) = 


F(z) 


r(l-a) 

One might consider that the right hand side arises from transforms having 

1-a 

values z and F(z). Now 
L(t /? ] = z _(1+/ * J r( 1+0 ) 

and if one chooses 0 so that -(1 + 0) = 1 -a, i.e. 0 = -(2 - a), then the 


integral for L[t ] diverges at the lower limit. Hence, as in Example 9.7.2, 
we consider 

z 1 


U(z) = 


r(l-a) a 
z 


• F ( z ) . 


Now L[ ■ 


1-a 


■T(a) 


(taking 0 = -(1 - a) above) and so, by the convolution 


theorem, 

T -l r r(a) 


[—-——• F( z ) ] = f 

a J 


f (r ) 


0 (t-r ) 


1-a 


dr , 


-1 


where L denotes the inverse operator. Consequently, as in Example 9.7.2, 

, t 


u(t) = 


f (r ) 


r(a)r(l-a) dt 1-a 


dr 


Sin net 


f (r ) 


dt J 1-a 

0 (t-r ) 


dr . 


17. Let U , U be the Laplace transforms of u , u ; K , K , K , K be the 

x JL x X x O 1 

transforms of k^, k^, k^, k^; and F, G be the transforms of f, g. Then 

U ! = F + + K 2 U 2 , 

U 2 = ° + K 3 U 1 + K 4 U 4- 
and hence 

(1 - K 1 )U 1 + K 2 U 2 = F ' 

- Vl + (1 - V U 2 = G ‘ 

One now solves for U^, and then applies the inverse transformation. 
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Chapter 10 

Asymptotic Expansions 


| READ §10.1 | 


At the top of p.438 one really requires more than 
h. 


k+1 w 


^ 0 . 


k+1 


What one requires is that for each k there is a neighbourhood of z^ in which 


k+1 w 


is bounded away from 0. From the conditions stated, in some 


k+1 


neighbourhood of z^, 


w 


h, la, , 

i < — 


k+1 


and hence in this neighbourhood, 


h. 


k+1 w 


k+1 


> a 


k+1 


h. 


k+1 


k+1 


Self-Assessment Questions 

SAQ 1, 2, 3 are Exercise 10.1, #1, 3, 4. 


| READ § 10.2 | 


Though the results of this section may not be unexpected the details of the 
proofs can be rather tricky. 

Problems for §10.1, §10.2 

P.1, 2, 3 are Exercises §10.1, #2, 5, 6. 

P.4, is Exercise 10.2. 


| READ §10.3 | 


The details in this section are a bit overwhelming, but one should try to 
get hold of the ideas. I feel that being adept at the kind of calculations 
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that arise is a bit like being highly trained for some sporting event. Most 
of us for most of the time are a bit flabby, but with a bit of common-sense 
and effort we can attain our goals if we take our time, and when it is 
necessary we can get into tip-top shape with a bit of training. 


Theorem 10.3.1 In this theorem we have conditions on F(t) for t small and 
for t large. We consider for <5 > 0 small, 


00 5 00 

F(t)e dt = J F(t)e Z ^dt + J F(t)e Z ^dt 
0 0 5 


8 00 

and use the two conditions on F(t) in J and in J respectively. In the 

0 8 


00 


8 oo 


proof of the theorem this splitting of J 

0 


into J and J 
0 & 


is not explicit, 


but it is coped with by the result for fixed N stated at the top of p.447. 
This incorporates the two conditions of the hypotheses. However, the 
splitting should be kept in mind as it may be useful in other contexts. 


Example 10.3.1 It is perhaps a bit misleading to say near the bottom of 
p.447 that Watson's lemma is being applied to the integral above. This 
integral requires quite a bit of modification before Watson's lemma can be 
applied. 


p . 448 Near the bottom of the page one requires that a , s are bounded away 
from 1, not just that a , s ^ 1 for u > s > 0. It is clear, in fact, that 
the stronger result is true for u > e >0. 

pp.448, 449 The explanation given for the transformation relating ? , u is 
somewhat obscure. One might be led to such a transformation by the 
following considerations. It is usually simpler to have a particular point 
at the origin and so we might see the effect of setting s - 1 = u . We are 
aware that we are looking for two solutions of u = s - 1 - log s and so we 
2 

might put 2u = ? and hope that these solutions will be coped with by the 
index 2? or rather we would note that for small w, 

2 3 

W - log(1 + U ) = -- — + ... 

2 

and then the appropriateness of 2u = ? is clear. 
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Theorems 10.3.2, 10,3.3 One should aim at acquiring an understanding of 
these results in broad outline rather than spending a lot of time over the 

oo S oo 


details. 

theorems. 


Note that one considers J 

0 



in the proofs of these 


p.452 At the top of the page the reference to Theorem 10.3.2 is perhaps 
inadequate. One has also to take account of the use of integration by parts 
as presented at the top of p.450. 

p.453 The description of the steepest descent method is more or less 
correct, but the assertions made are not so obvious as the presentation 
might lead one to expect. The method requires one to look at the zeros of 
h’(f) and show that through such a zero, ? say, one can find an admissible 

path of integration on which Re[h(?)] decreases as quickly as possible from 

Re[h(fg)]. For large x > 0, |e | will be relatively small on such a 

path except near to ? and so one hopes the behaviour of the integral can be 

found from Watson's lemma. More generally one hopes this idea will work for 
Re z > 0 and appropriately restricted. The steepest descent method as 
presented here has modifications which extend its usefulness, but we shall 
not go into this topic. 


Consider a point where h'(?) has a simple zero. We shall not consider 
other possibilities. Assume, without loss of generality, that ? = 0, 

h(?^) =0. If ? = pe 1 ^, then near ? = 0, 


h(? ) 


2 3 

a 2 t + a 3 f + 

■ | 2 i ( 2(/>+a ) 

lajp e 


+ 


(a ] _ * 0) 

(a = arg a 2 ). 


If a. = a„ = ... =0, then 
3 4 


Re[h(<r ) ] 



cos ( 2<j> 


+ a ) 


and the steepest descent will be when 2 <f> + a = ±n , i.e. <j> = —-—, 

<j> = - -—Ip general, when p is small the steepest descent path will 
leave ? = 0 in these directions. 
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Consider next a point ? on the path of steepest descent. We can take ? = 0 
and we suppose that h'(0) * 0. Then, if ? = pe ^, near ? = 0, 
h(? ) = b Q + + t> 2 ? 2 + ... (b^ 0 ) 

= b 0 + |b 1 1 pe 1 ^^ 5 + ... U3 = arg b 1 ). 

If b 2 = b 3 = _ = 0, then 

Re[h(f)] = Re b Q + |b 1 1 p cos(«£+/}) 

and the path of steepest descent is when <t> + 0 = n , i.e. <j> = n - . Now 

Im[h(f )] = Im b Q + | b^ | p sin (<p+/3) 


and 


Im[h(pe 1 ^)] = |b | sin(<£+/?) 

dp 1 

= 0 (<f> = IT - /? ) . 

One sees in general that if the path of steepest descent leaves { = 0 in the 
direction then 


_d_ 

dp 


Im[h(p e 1 ^ ) ] 


p = 0 , <t>=<f> 


0 


0 . 


Hence, if ? = ?(s) is the path of steepest descent parameterised by arc 

length, then 


ds 


Im[h(<r (s) ) ] 


and so Im[h(? )] is 


= 0 , 


constant on the path. 


Example 10.3.3 The wording of this example is somewhat odd. It is strange 
that arg z is not considered till later. An example of such complexity as 
this seems rather inappropriate as an illustration of the steepest descent 
method. 

p.457 One should consider the method of stationary phrase with caution. 

Its formal application may give the 'right' answer, but it may not be easy 
to justify the result. The conditions considered here are rather special. 


Problems for sl0.3 

P5, 6 , 7, 8 , 9, 10, 11 are Exercises 10.3, #1, 2, 3, 4, 5, 6 , 7. 


| READ § 10.4 | 
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This is a difficult section to get to grips with. You should note that the 
arguments are somewhat heuristic and justifications are adumbrated rather 
than dealt with properly. 


pp.461, 462 It is not made clear what distinguishes normal solutions from 
subnormal solutions. It would seem that in the former ^ 0 and in the 

latter = 0 and that is all. Otherwise consider the forms of the 

solutions in the two cases? 

p. 462 You should note that 7 ^, 7 2 are not arbitrary constants. They depend 

on the solution u. In fact 

2a z -2(5 du 
7, = e z 3— 

1 dz 

“~~1 


- "2 * u(z i>- 


Also in the definition of K(z, 7 /), z^ should be z. 

p.463 The discussion preceding Example 10.4.1 is confused. It seems to 
indicate that one can treat the case 7=0 and supplement the result in 


— 2^ ^ 2 j() 

this special case by adding on the term - 7^ J e ? d? . However, it 

z 

does not seem this can be done except, perhaps, in particular examples. 
p. 466 In the integral equation for 7? (4 ) the initial conditions should be 
c 2 = 7i (0), c ± = n ’ (0)/A . 


It is not clear to me how the integral equation is obtained as a consequence 
of the results in Exercises 7.2. However, one can argue as follows in order 
to obtain the equation. Given the initial conditions there is a (unique) 
solution 7? (4 ) of 

& 2 v 2 

—- + A 7] = r U )»f • 

a$ 

Let 7i (£ ) denote this solution and consider the equation 


. y* — \ S / 'l f 

a* 

i.e. we consider this equation for w with r(( )t| = r(4) r /(0 as a given 
function of <. The complementary function, i.e. the solution of the 
homogeneous equation, is 

w(0 = A sin A 4 + B cos A( (A, B constants). 


If D 


ci 

— then a particular integral comes from 
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w(£ ) = 


2 2 
D + A 


UU )nii)] 


(: 


2iA V D - iA D + iA 

4 


r-)[rU )» lit)] 


= — [ f 

2 U J 0 


e iA ({ t) - e 1A (? 11 ] r (t )j? (t )dt 


J sinUU “ t)] r(t)i?(t)dt. 


Now 

w(0) = B, w'(0) = AA 

and if we impose the same initial conditions on w as on r? then w(£ ) - 
satisfies the homogeneous equation and w(0) - r/ (0) = 0, w'(0) - r? 1 (0) 
and so wU) = v U ) and hence v satisfies the integral equation. 


Problems for §10.4 

P12, 13, 14, 15 are Exercises 10.4, #1, 2, 3, 5. 

Solutions to Self-Assessment Questions 

1. It is enough to show that 

N-l 

£(Z) = l a n W n (Z) + ° (W N- 1 ) (Z * V' 
n=l 

since then one can work backwards down the sequence step by step. By 
hypotheses 

N-l 

f(z) = I a nV Z) + a N W N (Z) + ° ( V (Z * V 

n=l 

and, since (w ) is an asymptotic sequence, 
n 

W N = ° (W N-1> (Z *V- 

Hence 

N-l 

£(z) = l w n (z) + otw^p (z z Q ). 
n=l 

2. There is some number R > 0 such that 

00 , 

f(z) = l a (z - z ) (|z - z 0 l < R). 

k = 0 

R 

The series converges for |z - z | = - and hence there is a constant K 
that 


« ) 
= 0 


such 
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Ia k | (-) < K (k = 0/ If 2,...). 

R 

Hence for |z - Zq| = - and any positive integer N, 

oo oo „ k k-N-1 

I l a k ( Z - V I S lz - z Q \ I K( 5 ) .(j) 
k=N+l k=N+l 


N+l 


N+l oo , k 


< K.(|) |z - z | I (|) 

k = 0 


N+l 

= 2K.(-) |z - Z Q 


N+l 


Therefore, for 0 < |z — z^| < 

N k 

f(z) = l a (z - z ) + g (z> 

k = 0 


where 


lg N (z)I o N+1 

T - ^ Tn s 2k ’ ( S’ ' z - 2 o' 


z-z 


0 


-> 0 ( z + z Q ) . 

. . n ... 

Since the asymptotic sequence in this case is w n ( z ) = ( z “ z q) this proves 

the desired result. 

3. We have that 
N 

f (z) = l a w (z) + a W N+1 (Z) + o(w ) (z - «„), 
n=0 

and since 

a N + l W N + l + ° (W N + 1> = 0<W N + 1> 
the result follows. 


Solutions to Problems 
1. For each positive integer N, 

, n N 2N+2 

2 4 , , N 2N 1 + (-1 t 

1 - t + t -+ -1 t = ---hr- 


1 +t 


and so 


*oo Z t oo m 

J 5—2 at = J e~ zt [ i - t 2 + ... + (- 1 ) t ]dt + r n (z), 

0 1+t 2 0 


where 
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N+l 


V z > = <- x > I 


, 2N+2 -zt 
t . e 


dt. 


0 1+t 

The integrals all converge since 

Rez > |z| cos - a) = |z| sin a 


and 


-zt, , -Rez.t, -|z|sina.t 
e = e < e 


Using integration by parts, for any positive integer k, 


r k -zt,. 1. k -zt 

| t e dt = - -t e 
J „ z 

0 


k r k-1 -zt,, 
+ - J t e dt 

0 z 0 


k P , k-1 -zt,. 

- . J t e dt. 


Hence 


r k -zt,. k! r -zt, 

J" t e dt = T . | e dt 

0 z 0 


k! 


k+1 


Therefore 


J e Zt [l - t 2 + ... + (-1) N t 2N ] dt 
0 


1 2! 4! 

z 3 + 5 

z z 


+ (- 1 ) 


N ( 2N) ! 
2N+1 


Also 


|R n (z)| < 


oo 2N+2 

J -2 6 dt 

0 1+t 


OO I I 

- 2N+2 - z sina.t,, 

< j t e dt 

0 


(2N+2)! 


i i • ,2N+3’ 

( | zIsina ) 

Hence the result follows with 


V z) = 15TT- 

Z 

2. First of all note that in R, 

Rez > |z| sin a > 0, |z + t| > Re(z + t) > \z\ sin a (t > 0) 
and so the integral is well defined. Integrating by parts 
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z -t 
e 


z + t 


dt = 


-t 


z+t 


00 00 


-t 


- J 


dt 


i - [- -1_ 

(z+t) 2 


1 1 „ r 

= -— + 2 f 

z 2 J 

z 


0 0 (z+t) 

-t 


-t 


- 2 J 3 

0 0 (z+t) 


dt ] 


-t 


dt 


0 (z+t) 

and so on. Hence for N a positive integer, 

z e”*' 1 1' 2» (~1) N ^.N 1 

J ?vr<3t = r- ;! + + 4- -.. + 1 ~ ~~~~ + R„(z), 


. z+t ’ z 23 *’’ N+l 

0 z z z 


N 


where 


-t 


V z, = (-1, .Nl J ——J 

0 (z+t) 


dt. 


Also 


I r n < z >I s ni j 


-t 


0 |z+t| 

N! 


N+2 


dt 


i | N+2 - A 

(|z|sinft) 0 

N! 


J e t dt 


, . .N+2 , , N+2 

( sm<*) I z I 

and so the result follows. 

3. Integrating by parts 


00 -t -t 

j V dt - - V 

z 


oo oo — t 

'J =-2« 

z z t 


-z 


-t 


- J 


dt, 


z t 


where we consider the integral along a part of a ray from the origin in the 
complex t-plane. On the path of integration |t| > |z| = r and 

. -t, -111COS0 i6> 

e = e , z = re , and so 


-t 


dt 


z t 


1 r -p cose. 

< -j J e dp 

r r 


1 e 


-p cos<9 


COS0 


sina 


r 

-z 


‘T5T 
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and hence the result follows. 

4. We interpret 

f(z) - g(z) (z € R, z -> z Q ) 


as 


k JL -► l (z € R. z -* z.) . 
g(z) U 

In this case 
f ( z ) n 

g (z) 0 

and so [f(z)] n ~ [g(z)] n . 


Let R = (£, the complex plane, z^ = oo and f(z) 
f(z) ~ g(z) , but 

. e «(*)-g(*> . e * + 1(z -> c) 

e 9<z) 

f( z ) g(z) 

and so e is not asymptotic to e 


Z 2 + 


g (z) = z 


Then 


Let R = {z: |z| > l}, 


z =oo and f(z) = 1 + g(z) = 1 + 
0 z 


2 * 


Then 


f(z) ~ g(z), but 

log f(z) _ log (1 + 1/z) _ 1/z - l/2z +... ^ ^ ^ 
log g(z) - lQg (1 + 1/z 2j ' 1/z 2 _ i /2z 4 +.. 

= Z [ 1 + 0(1) ] (Z -► 00) 


{ 1 ( Z -> oo) 

and so log f(z) is not asymptotic to log g(z). 

5. This is rather a tricky problem whose solution we shall only outline. 


Replacing u by t on p.448 et seq we consider F(t) 


d<7 

dt 


ds 
dt * 


Then F(? ) is 


analytic in {Re ? > 0} and behaves 'nicely' in {|arg H < 2 8 ^ for any 8 


with 0 < <5 < n ~. Consider z = x > 0 and J e F(? ) df round the contour 

consisting of the segment of the real f-axis from 0 to R, the arc of 

{| ? | = r} from ? = R to ? = Re 1/? with 0 = - - |) and the segment of the 

ray arg ? = 0 from Re 1/? to the origin. Letting R -> oo the integral over the 
arc tends to 0 and so 
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oo i 3 . , 00 

J e” Xte F(te 1/? )e 1/? dt = J e~ Xt F(t) dt. 

0 0 

Replacing x by z with 0 < arg z < n - a and noting that 

,n a . n a 


,ir ol i(i . 

~(2 " 2> - arg [ze ^ - * 


( 


a _ tr 

V ~ 2 


2 — z 

we see that when multiplied by z e the integral on the left with x 
replaced by z continues the r-function into the sector {0 < arg z < n 


a }. 


For small t we find that F(te 1 ^) has the same kind of expansion as F(t) 

except that t is now replaced by te 1 ^. If we now use Theorem 10.3.1 with z 

replaced by ze 1/? then we obtain the desired result. The sector 
{-(ff-ft) < arg z < 0} is dealt with in a similar fashion. 

6. Since n! = r(n + 1) = n r(n) from the asymptotic formula for r(z) on 
p.449 we obtain 


. -n n-1/2 
n ! ~ n.\J2n . e n 


. n -n 
= v2n-n.n .e 


7. This is the integral given on p.447 since 

1-s 1-s+logs 
se = e 

if we make the correction to the sign of log t. We proceed in a simple way 
as in the proof of Example 10.3.1. Putting t - 1 = u we find that 


r e z(l-t+logt) dt = j. eZ (-u+log(l+u)] du> 
0 -1 


We now use the monotonicity of -u + log (1 + u) noted in Example 10.3.1 and 

2 


u 

-u + log (1 + u) = -u + u~ — + 


(Iu| < 1). 


Hence we obtain as z ^ « with |arg z| < - - a, 

f e z[-u+log(l+u)] du _ 2 j“ e '^ 2 /2 du 
-1 0 


2 -t 

= 2.J- f e dt 
z J 


= J 


2 n 


To evaluate the integral we can take z = x > 0 at first and then use 

2 

analytic continuation. If t = s, then 
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r ^ ,, 1 f —1/2 — s 1 1 ■\/;r 

Je dt=-Js / eds=- r(-) = 

0 0 ^ 2 


2 * 


8. Assume |arg z| < - - a. One shows as in P5 that 


J t Z e Z ^dt = J 

0 0 


ooe 


-ie 


Z — ZC 

$ e df (9 = arg z) 


where the second integral is along the ray arg ? = -e and t Z = e Zl ° 9? with 
log t analytic in the right half ? -plane and log <r is real for <r real and 
positive. Hence 

OO OO II 

r z -zt r -10 z - z ho -i0 

J t e dt = j (pe ) e ' .e dp 

0 0 


-i(z+l)0 r 

e J o 

0 


V |z| 'dp. 


Put |z|p = u and then 


J t e 
0 


z -zt,. -i ( z+1 )9 
dt = e 


r . u v -u du 

I ( FT> ' e 


1 r 2 

—— u .e du 

z+1 J „ 
z 0 

T(z+1) 

z+1 

z 

r(z) 


Hence the asymptotic expansion of the integral follows from that of P(z). 

9. Consider first of all 

oo 2 oo iu 

J e 1 dt = - J du. 

0 2 0 ^ 

The latter integral clearly converges at the lower limit and integrating 

oo iu 
r e 

J — 77 T du by parts shows that it also converges at the upper limit. 

1 

Consider next 

J e lf2 * 

round the contour consisting of the segment [O, R] of the real axis, the arc 
(Re 1 ^: 0 < <p < and the segment [Re i?r//4 , 0] of the ray {arg <• = -}. The 
integral over the arc is 
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7T i/ 4 2 

/ iR (cos2<?!> + isin20 ) . i<f> , 

e . iRe a<f >. 

0 

2 4 (b t r 

Now sin 2<j> > -,2<p = — for 0 < <p < - and so the modulus of the integral 

n n 4 

above is bounded by 

7t/4 n 2 • ff/4 .2 J . 

c — R Sin2</> r — 4R (p/n 

Rj e d <j> < Rj e d <j> 

0 0 

' n - R / 

= 4R (1 - 6 > 


0 (R ■> oo) 


and so 


°o .2 , . 00 2 

r it iff/4 r -u 

J e dt = e I e du 

0 0 


= e 


i?r /4 vV 


Therefore 


00 , 2 00 , 2 

C(x) + i S(x) = J dt - J e lt dt 

0 x 


y/n . 1 , i. r " i t~ 

= “2 ( 72 + 7?' - I 6 dt ' 


oo i 2 2 

In J e dt put t = u and then 
x 

00 2 oo iu 

I e ' dt = 2 I , -172 dU - 

X 2 u 

X 

Integrating by parts we obtain 

iu 


oo iu 

f du 16 


X 


2 u 


1/2 


i 1/2 
u 


oo oo iu 

1 r e 

, + 2l J “V2 dU 

2 2 u / 

X X 


. 2 

IX , iu 

. e 1 r e 

1 x + 2i . 3/2 


i .u 


00 _ 00 iu 

3 r e , , 

, + 2l J , “572 dul 

2 2 u ' 

X X 


and so on. After the above steps 

J 


iu 


2 2 

oo iu ix , ix „ oo 

e . e 1 e 3 r 

0 1/2 x 2 3 2 J „ -I—- 

2 u x 2 l 5/2 

x xu 


du. 


Hence 
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00 


C(x) = \ ~ Re[ 5 J 


1U 


2 J 1/2 
2 u 
x 


du] 


2 2 

1 Jr 1 Sin X 1 COS X 

= 2 ^ " 2 1 - x + 2 -3 


cos u 


2 u 


5/2 


du], 


with a corresponding result for S(x). 


For any positive integer n, 

oo iu oo 


du 


'I , ~5Ti72 du| £ J , TTI72 

2 u 2 u 

x x 


2n-l 2n-l 
x 

Consequently the result of repeated integration by parts together with the 
inequality above shows that the series of the problem are asymptotic series 
10. Note from p.316 that 

. , 1 / -in<? -ysinfl _ 

J (iy) = J e * 69. 

n 2 7T J 

—n 

When y is large and positive, then one expects that the integral will be 

7T 

dominated by the contribution over a small interval about - -. Put 

e = - - + 6 and consider 
2 

1 r 3n//2 -in(-»r /2+<t> )+ycos</> 

d ( 1 y ) = n J e d< ^ 

n 2;r -n /2 


. n 3?r/2 


= 57 J 


-in^+ycos<^> , 
e d <p . 


-n / 2 


2 4 

Since cos <j> = 1 ~ Yl + JT • • • ' given e > 0 for all small 8 > 0 


1 - ( 1 + e ) — < COS <j> < 1 - (1 - e ) — ( k | < O • 

3tt 

For - - < <p < -8 and for 6 < <j> < assuming 8 is small enough, 

2 2 


.20 _ ,26 ^ 6 
cos 0 = 1 — 2 sin - < 1 - 2 sm - < 1 - 2 (-) 

Z /* Tl 


Also, for \<t> | < 8 , 


-in<£ 


- 11 =2 | sin < n5 . 


The proof of the asymptotic formula requires a number of steps. 
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(i) |<j"' + J 3 ' 72 ) e- in ^ +yCOS0 d^ | < 2,.e y[1 - 2(SA)21 

—n /2 5 

= 2,e' 2yU/ ' ,2 .e y . 

p -in0+ycos<£ r ycos<£_ 4 r -in<£ ycos<£ 

( 11 ) J e d<£ = J e dcf> + j (e - 1) .e d</> 

-8 -6 -8 


and 

i; (e _in ^ - l)e yCOS ^«| < n6 / 
-8 S 


e yC0 ^d^. 


(iii) Consider next 


y r „-y(i+e )4> 2 / 2 d < / ^ycos<^. 


e- J e 
-8 


'd <£ < J e J 


d* < e y / e -y(l->* 2 /2 


d<£. 


(iv) / e' y( 1 +£) ^ 2 / 2 d^ = _f e ”Y (1+£ >^ 2 / 2 d^ + A, 

—8 —oo 

where 


A = 2 J e- y(1+£ >* /2 d* 
8 


< 2 J e * y(1+£)S A/2 d^ 

6 


4 . e -y(l+ e )* /2 


y(1+e ) .8 


(v) j“ e _y(1+£)02/2 d« = 2 j“ e - y(1+£) * 2/2 d* 
—00 0 


= V- 


2tt 


y(l+e ) ' 

as in the solution to P7. 


I will not put everything down in detail, but indicate from the above how 

one arrives at the desired answer. Considering the integral over [- -, —J 

for J n (iy) one sees from (i) that for y > 0 large we can restrict attention 

to [-8, 5] for any small 8 > 0. (ii) shows that if 8 is small then the 

effect of replacing e in< ^ in the integral over [-5 , 8 ] by 1 is 'negligible 1 
(iii), (iv) and (v) shows that when 6 is small and y is large the effect of 
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8 oo 2 

replacing J d<j> by e Y J e ^ 2 

—8 —oo 

of the above is that as y ■> oo, 


d<^ is 'negligible'. The consequence 


e' 7 * / 2 a* 


. n y oo 
i e 


.n y 
1 e 

\/2n y 

11. We assume that in J e t Iy dt the integration is along the ray in the 

z 

complex t-plane given by arg t = arg z. Integrating by parts one obtains 


-t. -v 


-t. -v 


; e t dt = - e t - v J e . t 


-t -U-1 


dt 


-Z -v 


= e . z - [ -v e . t 


-t -i/-l 


-t . -v-2 


- i/ (i/ + 1) J e .t dt] 


e [z v + z/z v + v {v + 1) J e ^.t V 2 dt. 

z 


[Note that a^ = 1 and not 0 as stated in Dettman]. 


After an appropriate number of integrations by parts we have to deal with a 
term involving 

00 

r -t -i/-n 
Jet dt 

z 

and 


|J e fc t V n dt| < |e Z | J t V n dt 


.-i/-n+l i -z 
z . e 


i^+n-1 

Therefore the result follows. 

r 1/2 J p 


12. From w = ue 


we obtain 


w' = (u' - |pu)e 1//2 5 P , 


w" = [(u" - |p'u - -pu') - | p(u' - |pu)]e 1//2 J P . 

Using these results for w, w', w" in w" + pw' + qw = 0 gives the equation 
for u. 


Ill 



From the solutions to P3 and P7 in Chapter 7 we see that the relationships 
given for behaviour at <x> do in fact hold. 


13. From z = and w(z) = - ■ we obtain 


W'(z) = 

( ? 

U' (f ) U(J ) 


2? 


4 ' 


_ r u "(?) 3 u'(0 U'(?) 4u(<r ) dz 

[ } ~ ~ 7 ^~ 2 * 4“ 4~ + 5~ ( dT } 

2? <r <r <r 


-l 


U" (<r ) 5 u ' (? ) 2u(f ) 

4 4* 5 6 ’ 


4f 


f 


Hence 


u"(t ) - 5 .itiLL + [-| + 4f 2 Q(f )] u(f ) = 0, 

* ( 2 


where 


q q 

Qti ) = q(? (q x ^ 0 ) . 


It follows that 


8 „ 2 n/ x „ 8+4q 2 
~2 + 4? Q(f ) = 4q x + -— + . . . . 

t $ 

Applying the argument on p.461 to the equation for u(? ) we see that in the 
present case a * 0, 0 = 0 . This follows from the relations near the bottom 
of p.461 with q^, q^ there replaced by 4q^ ^ 0 and 0. Hence we get a 

'solution' 


u(f ) = e l c f 
k = 0 

and in terms of z and w(z) we find that 

n ./ 7 00 -U 

\/zw(z)=e Y c ] < f 
k=0 


and so 


w(z) 


, w 

a v Z r 

e x 

k = 0 


-l/2-k/2 

z 


[The answer given by Dettman seems to be wrong.] 

14. The argument required to establish this result is very similar to that 
used on pp.343, 344 to deal with the Fredholm integral equation. (Note, 
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incidentally, that in the formula in the middle of p.465 that verifies that 

2 

v(z) is a solution the 2 in the first expression should be »/ and the 1 in 
the second expression should be omitted.) 


Let 


H[f] = H[f](z) = J K(z, i,).- - dr?. 

z v 

In this notation, with v as specified on p.464, 

K 

f x = 1 + H[f], 

f 2 = 1 + H[f 2 ] = 1 + H[1 + H[f]] 

= 1 + H[l] + H 2 [f] 

= V x + H 2 [f], 

f 3 = 1 + H[ V]l + H 2 [f] 

= v 2 + H 3 [f], 


and so on. Assuming f is bounded it follows as on p.464 that 


| H[ f,| S " | H 2 i f ] | < 


and so 


|H n [f]| < i- * 0 (n oo) 

n! | , n 

I z | 

for any fixed z * 0. Therefore, for any such z, 

f n+ 1 ( 2 ) “ v r (z) = H n+1 [f](z) ->0 (n -► oo) . 

Consequently the sequences (f^) and (v^) converge to the same function. 


Uniqueness follows by considering two possible solutions v(z), V(z). 

start the iteration process with f(z) = V(z) then by the above 

f -*■ v (n oo) . But in this case f = V for all n and so V = v. 
n n 

15. We have that 


J (x) 


n=0 


, v n, „ 2n+v 
•1) (x/2) 

n !r (i/+n+i) 


_±_ ( *j 

r(u+i) * 


l 


/ , Ml 2n 
(-1) (x/2) 


n=l 


n! (i^ + l) (v+2) . .. (t/+n) 


If we 
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and 


00 


I l 

n= 


1 


(- 1 ) 


nl( m+1) 


(x/2 ) 2n 
. . . (t>+n) 


< 


1 

V + 1 


OO 


l 


(|x|/2) n 

n! 


= J_ (e ( UI/2) 2 


u. 


Hence the result follows. 
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